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1.0  INTRODUCTION 

Military  training  operations  involving  high-speed,  low-level  flights  are  routinely 
carried  out  by  all  the  flight  operation  commands:  SAC,  TAC,  MAC,  ATC,  ANG,  and  Air 
Force  Reserve  as  well  as  by  the  AFSC  for  various  research  and  development  programs. 
Aircraft  types  involved  in  such  operations  include  current  technology  single-  or  multi- 
engine  aircraft  operating  at  high  subsonic  speeds.  Operations  of  this  type  are  conducted 
on  specially  designated  Military  Training  Routes  (MTRs),  which  are  continually  changed 
for  a  variety  of  reasons.  Each  new  route  requires  an  environmental  assessment  to 
determine  the  community  noise  impact. 

The  computer  program  ROUTEMAP  calculates  the  noise  levels  on  the  ground  along  a 
route  corridor.  Program  ROUTEMAP  requires  the  Air  Force  planner  to  specify  the  nature 
of  the  flight  activity  for  each  MTR  in  question.  With  this  information,  ROUTEMAP  will: 

•  Compute  the  noise  levels  in  decibels, 

•  Provide  an  estimate  of  the  probability  of  high  annoyance,  and 

•  Rank  the  aircraft  in  terms  of  highest  noise  contributors  at  a  designated  point 
across  the  route. 

ROUTEMAP  is  a  menu-driven  program  that  runs  on  any  IBM  PC  or  PC-compatible 
computer  system.  ROUTEMAP  requires  MS  DOS  Version  2.0  or  later,  with  at  least  one 
megabyte  of  available  disk  space,  at  least  6^>0K  (kilobytes)  of  random  access  memory,  and 
8087/80287  math  coprocessor. 

Section  2.0  of  this  report  contains  a  description  of  the  model  formulation  and  the 
underlying  accumptions.  Section  3.0  is  a  description  of  the  model  computer  program. 
Section  is  the  user's  guide  and  gives  complete  instructions  for  using  the  program. 
Section  5.0  is  the  programmer's  guide.  It  contains  descriptions  of  all  subroutines,  data 
configurations,  and  programming  conventions.  Section  6.0  is  the  maintenance  guide  that 
describes  computer  code  modifications  required  if  the  format  of  either  NOISEFILE  or 
OMEGA  10  were  to  be  altered. 
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2.0  MODEL  FORMULATION 

Program  ROUTEMAP's  formulation  treats  operations  under  visual  and  instrument 
flying  rules  (VFR  and  IFR)  as  two  separate  models.  Operations  by  the  Strategic  Air 
Command  (SAC)  and  the  Tactical  Air  Command  (TAC)  are  indicative  of  the  range  of 
differences  seen  in  MTRs  within  the  United  States. 

SAC  low-altitude,  high-speed  training  missions  are  always  conducted  under  instru¬ 
ment  flying  rules  (IFR).  These  missions  provide  terrain  avoidance  training  during  flights 
to  bomb  scoring  ranges.  Operations  follow  the  route  centerline  with  a  high  degree  of 
precision  and  have  terrain-following  stage  lengths  of  between  150  and  250  nautical  miles 
while  the  minimum  altitude  above  ground  level  is  300  to  600  feet.  SAC  operates  several 
dozen  routes;  these  routes  can  range  in  length  up  to  1,000  miles. 

TAC  conducts  missions  under  visual  flight  rules  (VFR).  The  objective  of  these 
missions  is  to  practice  low-altitude  point-to-point  navigation  within  a  designated  MTR. 
TAC  low-altitude  operations  have  stage  lengths  of  200  to  300  nautical  miles  during  which 
the  minimum  altitude  above  the  ground  varies  between  100  and  500  feet.  These  flights 
tend  to  be  within  a  few  hundred  miles  of  home  base  and  are  controlled  by  the  local 
command.  TAC  has  several  hundred  routes  between  100  and  300  miles  long. 

Program  ROUTEMAP  models  both  VFR  and  IFR  routes  by  treating  an  individual 
flight  track  as  a  point  source  moving  along  a  line.  Since  flight  tracks  along  an  MTR 
seldom  follow  the  exact  same  trajectory,  an  MTR  can  be  viewed  as  a  distribution  of 
moving  point  sources  with  the  resultant  cumulative  noise  level  at  the  ground  being 
of  interest. 

Measurements  on  several  MTRs  (References  1  and  2)  have  established  that  a 
Gaussian  distribution  in  the  horizontal  plane  is  the  distribution  that  best  describes  the 
spatial  activity  along  an  MTR.  The  impact  of  flight  track  dispersion  in  the  vertical  plane 
on  sound  exposure  level  has  a  minimal  and  often  negligible  effect  compared  with 
dispersion  in  the  horizontal  plane.  For  purposes  of  the  present  model,  therefore,  vertical 
dispersion  is  not  considered;  thus  the  aircraft  tracks  are  distributed  laterally  at  a  constant 
altitude  above  the  ground. 

TAC  may  have  several  routes  in  the  same  geographical  area.  These  MTRs 
conceivably  could  overlap  or  follow  some  predominant  terrain  feature  that  is  separated  by 
several  miles.  To  model  the  various  geometrical  combinations  of  MTRs,  ROUTEMAP 
defines  a  working  space  called  a  corridor.  This  corridor  has  a  maximum  width  of  20  miles 
and  imposes  the  restriction  that  each  MTR  must  be  aligned  to  each  other  by  an  angle 
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not  to  exceed  15  degrees.  ROUTEMAP  will  model  any  combination  of  aircraft  type, 
flight  operating  parameters,  and  location  of  the  MTR  relative  to  the  corridor  centerline. 

RO'JTcMAP  will  only  permit  an  MTR  down  the  center  of  the  corridor  when 
modeling  an  IFR  route.  Limiting  the  operations  in  this  way  simulates  typical  SAC 
low-level  routes. 

The  noise  environment  from  MTRs  is  unique  in  several  respects.  Events  are  highly 
sporadic,  ranging  from  a  maximum  of  five  to  ten  per  day  to  a  minimum  of  a  few  (less 
than  ten)  operations  per  week  or  two.  This  differs  from  community  noise  exposure 
scenarios  where  operations  tend  to  be  continuous  or  somewhat  regular.  Individual  events 
are  also  different  from  typical  community  noise  sources:  the  combination  of  low  altitude 
and  high  air  speeds  results  in  noise  signatures  with  high  levels  and  short  durations  which 
can  result  in  very  rapid  onset  that  may  produce  a  startle  response. 

The  primary  noise  metric  used  in  this  model  is  the  onset  rate-adjusted  monthly  day- 
night  average  A- weighted  sound  level,  Lj^rnr-  ^^is  cumulative  noise  metric  was 
developed  in  a  study  that  uses  the  best  available  existing  knowledge  (Reference  3).  It  is 
based  on  an  integration  period  equal  to  one  calendar  month  with  the  highest  number  of 
monthly  operations.  Future  studies  may  result  in  a  change  in  the  integration  period,  but 
this  can  always  be  normalized  to  equivalent  monthly  operations.  For  this  reason,  program 
ROUTEMAP  requires  operations  on  an  MTR  to  be  quantified  in  terms  of  monthly 
statistical  averages. 

Two  other  noise  metrics  computed  by  ROUTEMAP  are  equivalent  continuous  sound 
level,  Leq,  and  day-night  sound  level,  Ljn.  The  L^q  is  the  level  of  a  steady  sound  which,  in 
a  stated  time  period  and  at  a  stated  location,  has  the  same  A-weighted  sound  energy  as  the 
time-varying  sound.  The  is  the  equivalent  continuous  sound  level,  in  decibels,  for  a 
24-hour  period  from  midnight  to  midnight,  obtained  after  the  addition  of  10  dB  to  sound 
levels  from  midnight  to  7:00  a.m.  and  from  10:00  p.m.  to  midnight.  Note:  ROUTEMAP 
computes  Lgq  and  Ljp  based  on  an  integration  period  of  a  calendar  month. 

The  metric  most  widely  used  for  noise-compatible  land-use  planning  is  L^jp.  Air 
Force  planning  policy  includes  interpretation  of  Ljp  in  terms  of  compatible  land  use. 
This  is  based  on  relationships  between  and  the  probability  of  highly  annoying  the 
population.^  In  Reference  3  it  is  recommended  that  L^p^nr  va^’Jes  along  MTRs  be  applied 
to  the  same  interpretive  criteria  as  Ljp  values  in  other  circumstances.  Based  on  this 
procedure,  ROUTEMAP  also  computes  the  probability  of  high  annoyance  corresponding  to 
the  computed  L^jip^p^. 
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3.0  PROGRAM  DESCRIPTION 

3.1  Computer  Program  Capabilities  and  Features 

”  lie  capabilities  and  features  of  ROUTEMAP  are  as  follows: 

•  Calculates  the  and  across  a  flight  corridor.  Up  to  20  different 

aircraft  operational  conditions  can  be  specified  within  this  corridor. 

•  Calculations  adjust  for  the  air  absorption  and  the  ground  attenuation  using 
aigorithms  that  approximate  those  employed  in  NOISEMAP.^ 

•  Calculates  the  probability  of  highly  annoyed  people.^ 

•  Allows  for  multiple  MTRs  across  a  corridor  that  has  a  width  of  +50,000  feet. 
The  centerline  for  each  MTR  can  be  positioned  anywhere  within  the  corridor. 

•  Displays  the  results  as  a  table  and  plots  the  corridor  noise  levels  on  any  user- 
specified  plotting  device.  Graphical  output  utilizes  PLOT88,  a  standard  com¬ 
mercial  subroutine  library. 

•  Ranks  the  contributors  to  the  total  noise  exposure  at  any  position  along  the 
corridor. 

•  Uses  NOISEFILE  and  OMEGA  10  to  compute  the  sound  exposure  levels.  Both  the 
data  base  and  the  program  are  Air  Force  standards. 

ROUTEMAP's  limitation  of  20  different  flight  configurations  is  due  to  dimensioned 
size  arrays.  These  arrays  can  be  increased  by  changing  appropriate  dimensions  as 
discussed  in  Section  5.0.  The  NOISEFILE  data  base  can  be  increased  from  an  existing  243 
number  of  entries  to  300  entries  with  the  existing  dimen;ions. 

3.2  General  Description  of  Program  Organization 

ROUTEMAP's  organization  is  shown  in  Figure  1.  The  arrows  indicate  the  order  of 
access  to  the  various  routines  and  the  overall  flow  of  the  program.  To  run  ROUTEMAP 
requires  three  basic  steps:  define  the  operations  across  the  corridor,  calculate  the  sound 
pressure  levels,  and  display  the  results.  All  three  of  these  steps  are  controlled  by  the 
main  program  and  initiated  from  the  MAIN  MENU. 
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MAIN  PROGRAM 


Figure  I.  Program  Organization. 
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To  define  the  operations  across  the  corridor,  the  MTRs  must  first  be  drawn  on  a  map 
using  the  navigation  points.  The  user  must  determine  which  Air  Force  command  flies 
these  routes  and  the  kinds  of  operations  that  are  performed.  Next,  the  user  must  draw  a 
corridor  that  has  a  width  less  than  20  miles.  The  flight  tracks  must  be  inside  the  corridor 
and  be  aligned  to  each  other  by  an  angle  not  greater  than  15  degrees.  See  Section  4. 2.5. 3 
for  a  description  on  adding  multiple  corridor  noise  metrics  if  a  route  is  skewed  by  more 
than  15  degrees. 

The  user  is  now  ready  to  run  ROUTEMAP.  Pick  a  point  along  the  corridor  centerline 
for  which  the  noise  levels  at  the  ground  perpendicular  to  that  point  are  required. 
ROUTEMAP  will  step  through  a  series  of  window  routines  that  requires  data  input  to 
describe  the  operations  for  the  corridor  cross-section.  When  all  the  operations  have  been 
entered,  ROUTEMAP  returns  to  the  MAIN  MENU. 

Calculating  the  noise  levels  begins  with  ROUTEMAP  writing  a  file  to  disk  that  will 
serve  as  an  input  deck  to  OMEGA  10.  ROUTEMAP  executes  OMEGAIO  via  a  shell  routine. 
OMEGA  10  calculates  the  sound  exposure  levels  (SEL)  for  the  aircraft  specified  and  writes 
the  results  to  a  file  on  disk.  ROUTEMAP  reads  this  file  and  uses  these  SEL  values  when 
calculating  the  noise  levels  across  the  corridor.  At  the  completion  of  the  calculations,  a 
table  is  displayed  showing  the  original  input  conditions  and  a  second  table  displays  the 

^dnmr  route. 

To  display  the  results,  the  user  can  either  plot  the  data  or  rank  the  contributors  to 
the  total  noise  level  at  any  point  across  the  corridor.  The  plotting  routines  will  plot  the 
noise  levels  and  the  probability  of  high  annoyance.  The  devices  available  are  specified  in 
the  PLOT88  user  manual  and  includes  virtually  all  popular  printers  and  plotters. 
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4.0  USER'S  GUIDE 


This  section  describes  the  procedures  for  properly  installing  and  running  ROUTE- 
MAP.  Topics  addressed  in  this  section  include  a  description  on  the  type  of  data  an  Air 
Force  planner  needs  to  gather  prior  to  running  ROUTEMAP,  and  a  discussion  on  how  to 
interpret  and  present  results. 

4.1  Getting  Started 

ROUTEMAP  is  delivered  on  two  disks.  One  of  the  disks  contains  ROUTEMAP's 
source  code,  a  batch  file,  and  the  executable  version  of  OMEGA  10.  On  a  separate  disk  is 
a  copy  of  NOISEFILE*  renamed  as  NOISE. 

Before  installing  ROUTEMAP,  be  sure  that  the  host  machine  has  1  megabyte  of 
available  disk  space,  640  kilobytes  of  random  access  memory,  and  a  8087/80287  math 
coprocessor.  Also  resident  on  the  host  machine  must  appear  a  copy  of  Microsoft 
FORTRAN,  a  SPINDRIFT  library,  and  a  PLOT88  library.  These  are  commercially 
available  packages;  see  Section  6.3  for  manufacturers'  addresses. 

PLOT88  is  a  graphics  subroutine  library  patterned  after  industry  standard 
CALCOMP  style  calls.  If  a  different  graphics  package  with  similar  structure  is  available, 
an  experienced  programmer  should  be  able  to  easily  adapt  ROUTEMAP's  graphics 
calls  to  it. 

To  install  the  software  it  is  recommended  that  a  new  directory  be  opened, 
preferably  on  a  hard  disk.  This  is  accomplished  by  using  the  following  DOS  commands: 

MKDIR  (directory  name) 

CD  (directory  name) 

Example:  MKDIR  ROUTE 
CD  ROUTE 

Next,  copy  MAIN.FOR,  OPS.FOR,  CALC.FOR,  PLT.FOR,  RANKODR.FOR,  ADDER.FOR, 
SPINATTR.HDR,  ROUTE.BAT,  OMEGA10.EXE,  and  NOISE  from  the  two  disks  into  the 
newly  created  directory.  To  do  this,  put  the  disk  into  drive  A,  then  type  the  following 
DOS  command: 

COPY  A;*.*  C: 

*  These  are  special  versions  of  NOISEFILE  and  OMEGA  10,  configured  to  operate  with 
ROUTEMAP,  and  are  not  directly  compatible  with  operation  of  NOISEMAP. 
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Here  it  is  assumed  that  C  is  the  destination  drive.  After  typing  this  command,  type  i 

ROUTE  to  begin  the  compilation  and  linking  procedures.  The  new  file  that  will  appear  on 
your  directory  is; 

ROUTEMAP.EXE 

This  procedure  completes  the  installation  of  ROUTEMAP. 

4.2  Running  ROUTEMAP 

To  execute  the  program,  the  user  must  type  the  program  name  ROUTEMAP.  After 
pressing  RETURN  the  screen  will  appear  blank  and  the  program  will  begin  an  initializing 
procedure.  This  procedure  will  take  approximately  30  seconds  on  a  PC-class  machine  and 
about  10  seconds  on  an  AT-class  machine.  Following  this  procedure,  a  banner  is  displayed 
and  the  user  is  instructed  to  press  ANYKEY. 

4.2.1  Main  Menu 

ROUTEMAP  is  a  menu-driven  program  that  interacts  with  the  user.  The  program 
structure  is  centered  about  a  MAIN  MENU  (see  Figure  2)  that  has  seven  options.  Not  all 
of  these  options  are  available  to  the  user  at  any  given  time.  ROUTEMAP  highlights  those 
options  which  can  be  activated.  There  are  three  main  steps  to  running  program 
ROUTEMAP;  first,  the  operations  are  defined;  second,  the  noise  levels  are  calculated; 
and  finally,  the  results  are  displayed.  Each  of  these  options  can  be  activated  from  the 
MAIN  MENU.  A  description  of  each  MAIN  MENU  option  follows: 

1.  Define  Operations  On  Military  Training  Route  -  Through  a  series  of  window- 
driven  routines,  ROUTEMAP  queries  the  user  to  determine  for  each  aircraft  that 
uses  the  route: 

Whether  the  mission  is  by  instrument,  visual  flying  rules,  or  special  operations; 

-  The  aircraft  type; 

-  The  aircraft  operating  power; 

-  The  vehicle  velocity; 

-  The  altitude  above  ground  level;  and 

-  The  number  of  daytime  and  nighttime  monthly  sorties. 

With  this  information,  program  ROUTEMAP  returns  to  the  MAIN  MENU.  The 
menu  display  highlighted  items  now  include  Option  3. 


S 


Figure  2.  MAIN  MENU. 
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2.  Change  Default  Settings  -  Only  one  screen  appears  for  this  option.  The  default 
settings  are; 

-  Temperature, 

-  Relative  Humidity,  and 
Number  of  Days  in  a  Month. 

If  a  default  setting  is  unsatisfactory,  a  change  is  made  by  pressing  RETURN  and 
making  the  appropriate  update.  Note  changes  in  the  default  setting  will  only 
impact  the  noise  levels  if  Option  2  is  chosen  prior  to  pressing  Option  3. 

3.  Calculate  Noise  Level  -  This  option  performs  two  important  functions:  first, 

OMEGA  10  is  executed  using  the  flight  parameters  specified  in  Option  1  and  the 
results  from  this  program  are  read  into  ROUTEMAP;  second,  OMEGAIO  data  file 
is  then  used  to  calculate  the  noise  level  along  the  ground.  Results  from  this 
calculation  are  presented  in  a  table.  This  table  lists  the  noise  level  in  at 

distances  perpendicular  to  the  route. 

4.  Plot  Noise  Level  Metrics  (Use  Data  Stored  on  Disk)  -  ROUTEMAP  plots  on  any 
designated  device  both  the  noise  level  and  the  probability  of  high  annoyance. 
The  plot  abscissa  represents  the  corridor  width  and  is  bounded  by  the  program 
calculations.  Any  portion  of  the  corridor  can  be  plotted  provided  the  user 
remains  within  a  +50,000-foot  limit. 

5.  Rank  Order  Contributions  to  Total  Noise  Exposure  -  Selection  of  this  option 
gives  the  user  the  capability  to  rank  the  contributors  to  the  total  noise  exposure 
at  any  distance  perpendicular  to  the  route. 

6.  Decibel  Addition  —  Computes  the  sum  of  two  decibel  levels. 

7.  End  This  Program  -  Terminates  the  program  and  returns  the  user  to  DOS. 

4.2.2  Configuration  of  a  Route 
4.2.2. 1  Selecting  a  Training  Mission 

Figure  3  shows  the  first  menu  to  appear  after  selecting  Option  1  from  the  MAIN 
MENU.  The  purpose  for  this  menu  is  to  determine  the  standard  deviation  of  lateral 
dispersion  about  the  route  centerline.  Of  the  four  choices,  the  Air  Force  planner  will 
typically  pick  either  1,  2,  or  4.  Option  3  entitled  "Special  Operations"  is  only  used  if  the 
Air  Force  planner  intends  to  use  his  own  value  for  the  standard  deviation. 
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I  DSFINE  OPESiTIOHS  OH  HllITAHT  TBAIHIHG  BDDTE 

I 

I  Select  A  Traioinf  Missioo 

ID  Visual  f lying  Buies  | 

(21  Instruient  flriog  Buies  | 

(3)  Special  Operations 

(4)  Betriere  Preriously  Sared  Flight  Track  File 

Your  Selection  : 

============================ 


Figure  3.  Select  Training  Mission. 
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Option  1,  Visual  Flying  Rules,  are  operations  normally  conducted  by  the  Tactical  Air 
Command.  These  operations  tend  to  be  during  daylight  hours  where  navigation  is 
accomplished  visually.  Navigation  on  these  routes  is  point  to  point.  The  points  can  be 
official  segment  end  points,  prominent  landmarks  which  are  natural  to  follow,  or  "delta 
points"  defined  by  command  for  training  purposes.  Multiple  MTRs  can  exist  in  a  single 
route  corridor,  governed  by  a  combination  of  these.  Reference  2  contains  a  detailed 
discussion  of  this  type  of  operation. 

An  Option  1  selection  activates  the  screen  shown  in  Figure  4.  The  user  must  decide 
whether  the  route  has  a  single  dominant  MTR  or  is  composed  of  multiple  MTRs  in  the 
same  immediate  vicinity. 

Consider  as  an  example  two  VFR  routes  located  in  an  area  that  has  dominant  terrain 
features  separated  by  2  statute  miles.  The  Air  Force  planner  in  this  situation  would 
select  a  standard  deviation  of  1.25  miles  for  both  routes  and  establish  their  respective 
centerlines  2  miles  apart.  If,  however,  the  routes  were  located  above  cattle  pastures  with 
no  dominant  terrain  features,  the  planner  may  choose  to  represent  both  routes  as  a  single 
corridor  with  a  standard  deviation  of  2.5  miles. 

Strategic  Air  Command  (SAC)  conducts  training  missions  under  Instrument  Flying 
Rules  (Option  2).  There  are  two  general  types  of  missions  flown  by  SAC  under  IFR; 
missions  that  simulate  an  attack  on  a  single  site  with  several  runs  over  a  target  via  a 
racetrack  pattern,  and  missions  where  a  network  of  routes  pass  over  several  sites  to 
simulate  attacks  on  multiple  targets.  In  either  instance  these  missions  are  navigated  with 
a  high  degree  of  precision.  Typically  there  are  no  maneuvers  under  IFR  except  for 
navigation  corrections.  Operations  of  this  kind  have  a  lateral  dispersion  described  by 
Gaussian  distribution  with  a  standard  deviation  of  0.5  statute  mile.  Reference  1  contains 
a  detailed  description  of  this  type  of  operation. 

Retrieve  Previously  Saved  Flight  Track  File  (Option  4)  is  used  when  the  Air  Force 
planner  has  preserved  a  configuration  file  on  disk.  When  choosing  Option  4  the  user  is 
given  two  opportunities  to  correctly  enter  the  filename.  If  ROUTEMAP  is  unable  to  find 
the  file  specified  the  first  time,  a  message  "file  not  found"  appears  on  the  screen, 
followed  by  a  directory  listing.  Pressing  ANYKEY  causes  ROUTEMAP  to  ask  for  the 
desired  file.  The  user  can  force  the  directory  to  be  shown  by  deliberately  entering  an 
incorrect  name.  A  selection  of  Option  4  accompanied  by  the  successful  retrieval  of  data 
terminates  the  configuration  procedure.  The  user  is  then  able  to  proceed  from  the  MAIN 
MENU  in  calculating  the  noise  levels. 
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715DAL  FLIGHT  ROLES 


OPTION  111  -  Specific  sell  defioded  track! s):  standard 
deriatioo  is  1.25  statute  liles. 


or 


OPTION  12)  -  Nultiple  scattered  tracks:  standard  deviation 
is  2.5  statute  ailes  about  centerline. 


Your  Selection 


Figure  4.  Visual  Flight  Rules. 
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4. 2. 2. 2  Selecting  an  Aircraft 


The  "Available  Aircraft  From  NOISEFILE"  menu  (see  Figure  5)  is  activated  follow¬ 
ing  the  selection  of  a  training  mission.  NOISEFILE  is  a  data  base  containing  experimental 
noise  measurements  from  aircraft  operating  at  different  flight  configurations.  NOISE¬ 
FILE  contains  over  200  entries  of  which  there  are  in  excess  of  50  different  aircraft.  Each 
NOISEFILE  entry  contains  the  flight  parameters  under  which  the  experimental  measure¬ 
ments  were  recorded.  Key  parameters  are  shown  in  this  menu  to  guide  the  user  in  the 
appropriate  selection.  The  user  is  advised  to  carefully  examine  this  file  using  the 
scrolling  feature  and  make  a  selection  that  is  closest  to  the  flight  operation  under 
investigation.  OMEGA  10  will  use  the  selected  datum  as  a  reference  and  will  sccile  it  to 
the  operating  conditions  specified  in  the  next  menu. 

When  selecting  an  aircraft  the  user  may  scroll  through  the  file  using  either  the  UP 
and  DOWN  arrows  or  the  PGUP  and  PGDN  keys.  When  the  desired  ai  ;raft  and  operating 
configuration  has  been  identified,  the  user  records  his/her  selection  by  pressing  the 
ENTER  key. 

4. 2. 2.3  Input  Operation  Parameters 

Figure  6  shows  the  "Operation  Parameters"  menu.  This  screen  requires  the  monthly 
average  operating  condition  for  each  aircraft.  The  training  mission  and  aircraft  type 
appear  at  the  top  of  the  screen.  Depending  upon  the  mission  type,  the  screen  will  request 
input  for  five  or  six  entries.  "Centerline  Track  Offset",  shown  in  the  figure,  will  appear  if 
the  user  has  specified  either  VFR  training  mission  having  a  standard  deviation  of  1.25  or 
has  selected  "Special  Operations";  otherwise,  this  option  is  not  activated. 

ROUTEMAP  has  an  allowable  corridor  width  of  +50,000  feet.  A  centerline  track 
offset  must  be  within  these  limitations.  If  the  screen  does  not  request  user  input  for  an 
offset  value,  ROUTEMAP  assumes  the  MTRs  are  down  the  center  of  the  corridor. 

Realistic  entries  for  "Power"  and  "Speed"  are  essential  if  ROUTEMAP  is  to  yield 
correct  results.  Values  must  correspond  to  the  planned  operating  mode.  This  information 
can  be  obtained  from  the  route  designers  or  from  experienced  aircrews  who  fly  that 
aircraft  on  similar  operations.  The  power  and  speed  corresponding  to  the  datum  selected 
in  the  previous  menu  are  shown  in  parentheses.  THESE  ARE  NOT  DEFAULT  VALUES. 
Their  intent  is  to  show  the  order  of  magnitude  for  both  the  power  and  speed,  and  the  units 
in  which  the  data  are  to  be  entered. 


AVAILABLE 

AIRCRAFT  FROR 

NOISEFILE 

i 

POKER 

DRAG 

AIR 

AIRCRAFT 

DESCRIPTION 

POKER 

CONFIGDRATION 

SPEED 

1 

E-3A 

TAEEOFF  POKER 

1.83  EPR 

GEAR  DOKN,  50DEG  FLAPS 

250 

E-3A 

APPROACH  POKER 

1.45  EPR 

GEAR  DOKN.  5DDEG  FLAPS 

250 

E-3A 

IRTERREDIATE  POKER 

1.50  EPR 

SEAR  DOKN.  SPEED  BRAKE 

250 

E-3A 

TRAFFIC  PATTERN 

1.12  EPR 

NO  DRAG 

250 

EC-lOA 

TAEEOFF  POKER 

110.0  X  R1 

TAKEOFF  POKER 

230 

IC-lOA 

APPROACH  POKER 

79.0  X  R1 

APPROACH 

165 

Use  arrov  kers  to  scroll  desired  aircraft  into  highlighted  area, 
press  SlTDilS  shea  dooe. 


Figure  5.  NOISEFILE  Data  Base. 


?ISDil  FlieHT  HOLES 

oPERdTioR  numm 

Aircraft  Type  :S-3A  Coafiguratioa  :  1 

Poser  (  1.83  SPR  )  : 

Speed  (250  knots)  : 

Altitude  Above  Ground  Level  (feet)  : 

Centerline  Track  Offset  (feet)  : 

Ruiber  Of  Daytiie  Honthlr  Sorties  : 

Nuaber  Of  Righttiie  Ronthly  Sorties  : 

press  ESC  shen  all  done 


Figure  6.  Operation  Parameters 


When  entering  data,  the  user  inputs  a  value  and  presses  RETURN.  If  the  entry  is 
acceptable  ROUTEMAP  will  then  echo  the  number  and  advance  to  the  next  data  entry 
point.  As  a  check  that  all  entries  are  properly  recorded  by  ROUTEMAP,  the  user  should 
press  RETURN  successive  times  and  observe  the  cursor  movement  from  one  data  entry 
point  to  the  next.  If  all  is  well,  then  press  ESC  to  advance  to  the  next  screen. 

Two  actions  occur  when  the  user  presses  ESC:  first,  the  data  are  recorded  by 
ROUTEMAP;  second,  ROUTEMAP  flashes  the  screen  "Enter  Another  Configuration 
(Yes/No):".  A  "Yes"  entry  causes  ROUTEMAP  to  display  the  NOISEFILE  data  base  screen. 
A  "No"  entry  causes  ROUTEMAP  to  flash  the  screen  "Save  File  On  Disk  (Yes/No):".  If  a 
file  is  saved  on  disk  it  may  be  retrieved  at  a  later  time,  thus  circumventing  the  task  of 
reconstructing  the  corridor  data  base.  If  minor  changes  to  an  existing  file  are  desired, 
this  file  can  be  edited  with  a  text  editor  and  resubmitted  to  ROUTEMAP, 

Following  the  save  file  on  disk  screen  the  user  is  returned  to  the  MAIN  MENU.  This 
completes  the  configuration  procedure. 

4.2.3  Set  Default  Values 

Shown  in  Figure  7  is  the  only  screen  to  appear  after  pressing  Option  2  from  the 
MAIN  MENU,  The  default  setting  screen  contains  those  parameters  entered  into  the 
calculation  that  typically  do  not  change  from  one  operation  to  the  next.  To  change  a 
setting,  press  RETURN  in  succession  until  the  cursor  is  positioned  at  the  entry  to  be 
revised.  Now  enter  the  value  and  press  RETURN.  If  the  entry  is  accepted  by 
ROUTEMAP,  the  value  will  be  echoed.  When  satisfied  with  the  default  settings,  press 
ESC  to  return  to  the  MAIN  MENU. 

Note:  Setting  the  default  values  will  ONLY  impact  the  noise  levels  if  Option  2  is 
updated  prior  to  pressing  Option  3  from  the  MAIN  MENU. 

4.2.4  C2Uculating  the  Noise  Levels 

Figures  8  through  11  show  the  succession  of  screens  that  will  appear  after  pressing 
Option  3  from  the  MAIN  MENU. 

"RUNNING  OMEGAIO"  is  the  first  message  flashed  on  the  screen.  If  the  operation 
parameters  entered  during  the  configuration  process  are  realistic,  then  OMEGAIO  will 
successfully  terminate  without  an  error  message  being  displayed.  If  an  error  were  to 
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DEFiOLT  SETTINGS 

Teiperature  la  Degrees  Fahrenheit  :  59 
Relatire  Huiiditr  •'  70 
Nuiber  Of  Days  In  Hooth  :  30 

press  ESC  shea  all  done 
Figure  7.  Default  Settings. 
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min  mno 

Please  Hake  i  Selection 


(1)  Define  Operations  On  Hilitarr  Training  Route 


CkLCOLkTlNe  lOISS  LEVKLS 


(5)  Rank  Order  Contributors  To  Total  loise  Exposure. 

(6)  Decibel  Addition. 

(7)  End  This  Prograi. 

Tour  Selection  ;  I 


Figure  9.  Calculating  Noise  Levels. 


TBI  IRPDT  DATA 

Standard  Deviation  :  1.25  Statute  Biles 
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,0 
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5 

F-16 
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.0 
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85.00 

.528E44 

6 

F-5E 

15.0 

.0 
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-.528E+4 

7 

F-15 

20.0 

.0 
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-.528E44 

8 

F-14 
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.0 
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9 
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30.0 

.0 
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.000 

Figure  10.  Input  Data  Screen. 
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Figure  11.  Noise  Distribution. 
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occur,  ROUTEMAP  will  display  one  of  several  error  messages.  If  the  error  is  fatal, 
ROUTEMAP  will  terminate.  To  diagnose  the  nature  of  the  error,  examine  OMEGAIO.OUT 
file  and  check  to  see  if  a  file  IN  and  OUT  were  indeed  created. 

OMEGA  10  run-time  is  typiccilly  15  seconds  per  aircraft  type.  If  six  different  kinds 
of  aircraft  operated  on  a  given  corridor,  then  allow  about  a  minute  and  a  half  for 
OMEGA  10  to  calculate  the  sound  exposure  levels  (SEL). 

"CALCULATING  NOISE  LEVELS"  is  the  second  message  which  is  flashed  on  the 
screen.  At  this  time,  ROUTEMAP  is  calculating  the  actual  noise  level  on  the  ground  for  a 
corridor  having  a  width  of  100,000  feet.  The  run-time  for  this  procedure  ranges  from 
several  minutes  to  as  long  as  an  hour.  The  amount  of  time  required  is  dependent  upon  the 
number  of  routes,  their  respective  width,  and  the  number  of  configurations. 

At  the  completion  of  this  calculation,  a  screen  is  displayed  showing  the  original 
input  data  (see  Figure  10).  To  make  a  hard  copy  of  this  screen  for  your  records,  press 
PRINT  SCREEN.  To  continue  to  the  next  screen,  press  ANYKEY. 

Next,  ROUTEMAP  displays  the  noise  level  in  across  the  route  (see  Figure  11). 

As  before,  a  hard  copy  of  this  data  is  made  by  pressing  PRINT  SCREEN.  To  preserve  this 
data  on  disk  for  future  plotting,  answer  "Yes"  when  asked  to  save  the  file.  Otherwise,  the 
data  is  lost  at  the  termination  of  ROUTEMAP  or  if  a  new  corridor  is  defined. 

The  noise  distribution  table.  Figure  11,  contains  the  information  generally  needed 
for  preparation  of  an  Environmental  Assessment.  The  planner  is  usually  interested  in  the 
maximum  level  and  the  location  of  the  65  dB  contours.  In  the  example  shown, 

the  maximum  is  56  dB.  There  is  no  65  dB  region.  If  there  were  a  need  to  show 

Ldnmr  contours,  these  are  at  -14,000  feet  and  15,000  feet  relative  to  the  corridor 

centerline,  and  would  be  drawn  at  these  distances  and  parallel  to  the  centerline. 

4.2.5  Analysis  of  Results 
4.2.5. 1  Plotting 

Option  4,  "Plot  Noise  Level  Metrics",  activates  the  plotting  routines  from  the  MAIN 
MENU.  The  first  screen  that  appears  is  dependent  on  ROUTEMAP's  current  status.  If 
there  are  no  noise  levels  to  be  plotted  because  ROUTEMAP  has  not  made  the  calculations 
prior  to  the  selection  of  Option  4,  then  ROUTEMAP  assumes  that  the  file  to  be  plotted 
must  reside  on  disk.  In  this  situation  ROUTEMAP  queries  the  user  for  the  filename. 
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The  second  screen  to  appear  is  shown  in  Figure  12.  The  screen  is  divided  into  three 
sections:  General,  Noise  Level,  and  Annoyance.  In  the  General  section  the  user  inputs 
the  plot  title  and  the  range  of  distances  across  the  corridor  to  be  plotted.  This  subdivision 
is  entitled  "General"  because  these  parameters  will  appear  in  both  the  noise  level  and 
probability  of  high  annoyance  plots. 

Next,  the  user  must  choose  which  noise  metrics  are  to  be  plotted,  whether  symbols 
are  to  be  used,  and  if  all  of  the  curves  are  to  be  collapsed  onto  one  plot.  If  the  user  has 
selected  more  than  one  noise  metric  to  be  plotted,  it  is  recommended  that  symbols  are 
used  to  distinguish  the  curves. 

Finally,  the  user  indicates  whether  a  probability  of  high  annoyance  curve  is  to  be 
plotted  for  the  corridor. 

Once  the  user  has  completed  data  entry,  press  ESC  to  proceed  to  the  next  screen. 

The  third  screen  queries  the  user  to  select  a  plotting  device  (see  Figure  13).  Five 
devices  are  available  to  the  user  from  this  menu.  If  the  user  device  does  not  appear  on 
this  screen,  then  press  Option  6,  "Other".  The  user  must  now  enter  the  lOPORT  and 
MODEL  number.  The  user  should  consult  the  PLOT88  user  manual  to  ascertain  the 
appropriate  lOPORT  and  MODEL  number.  The  scale  factor  selection  is  done  by  trial  and 
error.  Note  that  the  present  devices  use  scale  factors  between  0.7  and  1.2. 

At  the  completion  of  the  third  screen,  the  plots  will  be  drawn  on  the  appropriate 
device  and  the  user  will  be  returned  to  the  MAIN  MENU.  The  user  has  the  option  to 
display  the  data  under  different  or  the  same  conditions  by  pressing  Option  4  from  the 
MAIN  MENU. 

The  plots  generated  by  this  option  can  be  useful  to  a  planner  in  developing  a  detailed 
understanding  of  the  route  noise  environment.  The  plots  are  also  suitable  for  inclusion  in 
an  Environmental  Assessment  if  the  presentation  would  be  enhanced  by  that  level 
of  detail. 

4. 2. 5. 2  Rank  Order  Contributors 

When  route  noise  levels  are  high  enough  to  cause  a  potential  problem,  it  is  useful  to 
know  the  relative  noise  contribution  for  each  type  of  operation.  This  information  can  form 
the  basis  of  noise  abatement  via  changes  to  route  operation. 
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PLOT  RESDITS 

GSMEBAL 

Enter  Title  : 

Enter  The  Distances  Across  HTB  To  Be  Plotted. 
Hiniiui  :-50630. 

Haziitis  :50D80. 

BOISE  LEVEL 

Do  Toy  Bant  To  Plot  Ldnir  (Tes/Bo)  : 

Do  Toy  Kant  To  Plot  Ldn  (Tes/Bo)  : 

Do  Toy  Bant  To  Plot  Leq  (Tes/Bo)  : 

Do  Toy  Bant  Sribols  Dsed  (Tes/Bo)  ; 

Do  Toy  Bast  To  Collapse  These  Curves 
Onto  One  Plot  (Tes/Bo)  : 

ABBOTABCE 

Do  Toy  Bast  To  Plot  The  Probahilitr  Of 
Bigh  Annorance  (Tes/Bo)  : 

press  ESC  shes  all  done 


Figure  12.  Plot  Menu. 


Figure  13.  Device  Menu. 


Option  5  from  the  MAIN  MENU  will  permit  the  user  to  rank  order  the  contributors 
to  the  total  noise  exposure  at  any  designated  position  across  the  corridor.  The  user 
specifies  the  distance  from  the  corridor  centerline  to  the  ground  location  point  of  interest 
to  be  ranked.  This  value  must  be  within  the  +50, 000-foot  limit.  Once  the  screen 
displaying  rank  contributions  is  displayed,  press  ANYKEY  to  return  to  the  MAIN  MENU. 


4. 2. 5. 3  Combination  of  Noise  From  Separate  Routes 

When  the  angle  between  two  routes  exceeds  15  degrees,  or  if  incompatible  opera¬ 
tions  (e.g.,  IFR  versus  VFR)  preclude  treating  them  as  a  single  corridor,  ROUTEMAP  must 
be  applied  separately  and  the  results  combined.  This  section  describes  the  method  used  to 
add  two  corridors  together  using  ROUTEMAP's  tabulated  results.  A  step-by-step 
description  follows: 


1.  Lay  out  the  routes  on  a  map.  Draw  the  corridors  and  respective  centerlines. 
Locate  on  the  map  the  ground  location  for  which  the  noise  level  is  to  be 
determined.  Call  this  point  the  receiver  position. 

2.  Draw  two  lines  on  the  map  which  are  perpendicular  to  the  corridor  centerline 
and  intersect  the  receiver  position. 

3.  Run  ROUTEMAP  for  each  corridor  to  obtain  the  noise  levels  at  points  perpen¬ 
dicular  to  each  corridor's  centerline. 


4.  PRINT  SCREEN  the  tabulated  results  from  ROUTEMAP. 


5.  Using  the  two  tables,  sum  the  noise  contribution  from  both  corridors.  This 
capability  is  provided  as  an  option  from  the  MAIN  MENU.  The  equation 
employed  is  as  follows; 


Noise  Level 
at 

Receiver  Position 


Number  of 
Corridors 

‘°8l0  X) 

i  =  1 


This  procedure  can  be  applied  to  the  analysis  of  a  particular  noise-sensitive 
receptor.  It  can  also  be  applied  to  a  number  of  points  around  two  corridors,  and  the 
results  used  to  plot  combined  noise  contours  for  use  in  planning  or  an  Environmentad 
Assessment. 
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3.0  PROGRAMMER'S  GUIDE 

5.1  Convention  and  Annotation 

ROUTEMAP  is  written  in  standard  Microsoft  FORTRAN  77.  All  real  variables  are 
single  precision,  REAL*4.  All  fixed  point  variables  are  INTEGER*4.  ROUTEMAP  makes 
extensive  usage  of  SPINDRIFT  and  PLOT88  library  routines,  commercial  packages  which 
are  PC  and  MS-DOS  specific.  A  PLOT88  Site  License  is  required  to  run  ROUTEMAP.  The 
program  source  and  executable  code  are  fully  transportable  to  any  IBM  PC-compatible 
machine.  Compilation  and  linking  of  source  code  to  produce  executable  code  (Section  6.3) 
requires  licensed  copies  of  Microsoft  FORTRAN  Version  4.01  or  later,  the  SPINDRIFT 
library,  and  the  PLOT88  library. 

Extensive  use  is  made  of  subroutines  in  order  to  provide  reasonably  well-structured 
code.  Variable  names  are  consistent  between  subroutines,  although  in  several  cases 
the  mnemonics  are  similar  but  not  identical.  Executable  statement  labels  begin  with  10 
in  each  program;  format  labels  begin  with  1000.  Text  strings  used  by  SPINDRIFT 
are  incorporated  into  the  code  and  window  identification  numbers  begin  with  1  in 
each  subroutine. 

5.2  Definition  of  Symbols 

Symbols  defined  herein  are  used  throughout  this  report  and  program  ROUTEMAP. 
A  sketch  of  the  MTR  and  the  use  of  the  symbol  notation  is  shown  in  Figure  14.  The 
symbols  listed  in  Table  1  are  a  complete  list  of  the  variables  used  in  ROUTEMAP. 
Dummy  variables  used  for  data  handling  and  integer  counters  are  not  included  in  this  list. 
Symbols  that  are  intended  to  be  arrays  will  be  listed  with  their  cirray  dimensions. 
Variables  I,  J,  K,  and  L  will  serve  as  running  indices.  Variables  I  and  J  are  usually  but  not 
always  used  to  represent; 

•  I  -  the  aircraft  index,  and 

•  3  -  the  offset  distance. 

I 
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50,000  Ft 


Figure  14.  Schematic  Illustrating  MTR  Symbol  Notations. 
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Table  1 


Symbol 

ACC(300) 

AIRSPDOOO) 

ALT(20) 

ATTEN 

BETA 

DIFF 

DRAG(300) 

FACT 

FINTEN(20, 147,3) 
FLDNM(201) 
FLDNMR(20l) 
FLEQ(201) 

FMAX 

FMEAN(20) 

FMIN 

GNDPOS 

ICODE(20) 

ICOL( ) 

lOPORT 


Variable  Definitions  Used  in  ROUTEMAP 


Description 

Aircraft  code  from  NOISEFILE. 

Aircraft  speed  from  NOISEFILE. 

Aircraft  altitude  in  feet  (see  Figure  14). 

Lateral  attenuation  in  dB. 

Elevation  angle. 

Distance  as  measured  from  center  of  the  Gaussian  distri¬ 
bution  to  point  along  distribution  that  is  presently  being 
evaluated  in  feet  (see  Figure  14). 

Aircraft  drag  configuration  description  from  NOISEFILE. 
Scale  factor  for  plotting  routine. 


A-weighted  mean  square  pressure. 

L^^  value  in  dB  for  nominal  monthly  integration  period. 
Ldnmr  value  in  dB  for  nominal  monthly  integration  period. 
Lg^  value  in  dB  for  nominal  monthly  integration  period. 
Maximum  X  axis  coordinate  selected  by  user. 


Flight  track  centerline  measured  from  reference  position 
centerline  in  feet  (see  Figure  14). 

Minimun  X  axis  coordinate  selected  by  user. 


Absolute  position  along  flight  track  with  respect  to  refer¬ 
ence  coordinate  in  feet  (see  Figure  14). 

An  integer  array  used  to  point  to  the  NOISEFILE  record 
that  the  user  has  picked  for  the  corridor  configuration. 

Cursor  column  position.  Closed  parentheses  are  used  to 
indicate  the  array  dimensions  vary  between  subroutines. 


INTEGER*2  variable  used  by  PLOT 88. 
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Table  1  (Continued) 


Symbol 

Description 

IROW( ) 

Cursor  row  position.  Closed  parentheses  are  used  to 
indicate  the  array  dimensions  vary  between  subroutines. 

ISWITCH 

Integer  variable  that  points  to  the  mission  type. 

MAX 

The  total  number  of  aircraft  records  stored  in  NOISE- 
FILE.  Not  to  be  confused  with  the  number  of  lines  in  file. 

MODEL 

INTEGER*2  variable  used  by  PLOT 88. 

MONTH(20) 

Number  of  days  in  the  month. 

NAME(300) 

Aircraft  name  from  NOISEFILE. 

NUMAIR 

Number  of  aircraft  type. 

NUMD(20) 

Number  of  daytime  monthly  sorties. 

NUMN(20) 

Number  of  nighttime  monthly  sorties. 

OFFSET 

Distance  from  receiver  location  to  point  along  Gaussian 
distribution  that  is  presently  being  evaluated  (see 
Figure  l^t). 

ONSET 

Onset  rate  in  dB  per  second.  The  Onset  rate  is  defined  by 
slope  of  the  A-weighted  time  history  between  points  5  dB 
above  the  ambient  and  5  dB  below  the  maximum  level. 

ONPEN 

Onset  rate  adjustment  in  dB,  which  is  applied  only  when 
the  maximum  A-weighted  sound  level  exceeds  the  ambi¬ 
ent  level  by  at  least  15  dB. 

OPCR(300) 

Operation  power  code  from  the  NOISEFILE  data  base. 

POWDES(300) 

Aircraft  power  description  from  NOISEFILE. 

POWSET(300) 

Aircraft  power  setting  from  NOISEFILE. 

POWUNT(300) 

Aircraft  power  setting  corresponding  units  from 
NOISEFILE. 

PROB 

Probability  of  aircraft  flying  at  specified  X  position  along 
military  route. 

RANK(20,201) 

L  .  values  for  each  individual  aircraft  at  500-foot 

increi^ients  across  corridor. 
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Table  1  (Continued) 


Symbol 

SCN(1024) 

SEL(20,25) 

SIGMA 

TITLE 

TOGGLE 

VEL(20) 

X 

XARRAY(203) 

XFINL 

XINIT 

YARRAY(203) 


Description 

Scratch  pad  used  to  construct  screens,  prior  to  their 
display. 

Sound  exposure  level  for  a  specific  aircraft,  power  set¬ 
ting,  and  airspeed,  in  dB. 

Standard  deviation  of  track  width  in  feet. 

Plot  title. 

Logical*!  variable  that  tells  the  plotting  program  whether 
the  user  has  called  the  main  plotting  menu  before. 

Aircraft  airspeed  in  feet  per  second. 

Distance  as  measured  from  reference  position  to  point 
along  Gaussian  distribution  that  is  presently  being 
evaluated  (see  Figure  14). 

Absolute  grid  coordinates  used  to  plot  X  axis. 

Upper  integral  limit  for  normal  distribution  calculation. 

Lower  integral  limit  for  normad  distribution  calculation. 

Absolute  grid  coordinates  used  to  plot  Y  axis. 
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5.3  Subroutines 


ROUTEMAP  subroutine  hierarchy  chart  is  shown  in  Figure  15.  The  main  program  is 
described  in  Section  3.2.  The  following  subsections  describe  each  subroutine,  presented  in 
the  order  they  appear  in  the  hierarchy  chart.  The  descriptions  here  and  in  Section  3.2  are 
intended  to  be  read  together  with  the  source  listing  in  Appendix  A  if  a  detail  understand¬ 
ing  of  the  code  is  desired. 

5.3.1  INITIAL 

This  subroutine  initializes  the  SPINDRIFT  library,  the  display,  and  the  keyboard. 
Subroutine  INITIAL  captures  the  present  screen  image  and  stores  it  in  array  INTSCN.  The 
keyboard  numberlock  is  cleareo.  The  video  mode  is  determined  and  set  to  80x25  color  if 
it  is  not  a  monochrome  display. 

5.3.2  OPEN 

This  routine  draws  the  opening  banner.  The  zoom  feature  is  accomplished  by 
repeatedly  drawing  windows  of  increased  size.  After  each  window  is  drawn,  a  two-tone 
harmonic  is  sounded  using  a  SPINDRIFT  library  routine.  A  call  made  to  subroutine 
PAUSE,  draws  in  the  lower  left-hand  corner  a  window  instructing  the  user  to 
press  ANYKEY. 

5.3.3  PAUSE 

This  routine  draws  a  window  in  the  lower  left-hand  corner  of  the  screen  with  a 
message  to  press  ANYKEY. 

5.3.4  MENUaVAL,FLAG) 

This  routine  draws  the  MAIN  MENU.  Logical  array  FLAG  is  passed  to  subroutine 
MENU  through  the  call  statement.  Array  FLAG  indicates  the  present  status  of 
ROUTEMAP.  If  FLAG  is  true,  then  an  option  appearing  in  the  MAIN  MENU  can  be 
activated.  Subroutine  MENU  uses  the  FLAG  array  to  set  the  brightness  attribute  when 
displaying  the  menu.  After  the  menu  is  displayed,  the  user  is  queried  for  a  selection.  If  a 
key  other  than  a  number  key  is  depressed  the  keystroke  is  ignored.  Following  a  successful 
entry,  IV AL  is  returned  to  the  main  program. 
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MAIN 
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Figure  15.  Subroutine  Hierarchy  Chart,  ROUTEMAP. 
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5.3.5  opsiaswrrcH) 

This  subroutine  determines  the  type  of  training  mission.  The  screen  entitled  "Define 
Operations  On  Military  Training  Route"  is  displayed  with  the  four  options.  The  computer 
waits  until  a  key  is  pressed.  If  the  key  pressed  is  not  a  number  key  then  the  keystroke  is 
ignored.  Out-of-range  entries  are  ignored.  If  the  user  requests  the  retrieval  of  a 
previously  saved  data  base  then  the  user  is  queried  for  the  filename.  A  "file  not  found" 
message  is  displayed  only  once  before  returning  the  user  to  the  previous  screen.  Integer 
variable  ISWITCH  is  passed  from  subroutine  OPSl  to  the  main  program. 

5.3.6  OPS2aSWITCH) 

This  subroutine  determines  the  operating  parameters  across  the  corridor.  The 
subroutine  is  divided  into  three  major  sections.  First,  subroutine  OPS2  queries  the  user 
for  the  standard  deviation.  Second,  the  aircraft  type  is  recorded.  Last,  the  operating 
parameters  for  the  designated  aircraft  are  determined.  The  program  iterates  between 
last  two  sections  until  all  aircraft  configurations  are  entered  or  twenty  configurations 
have  been  specified. 

Integer  variable  ISWITCH  passed  from  the  main  program  instructs  subroutine  OPS2 
which  training  mission  has  been  specified.  A  computed  GO  TO  statement  decides  which 
screen  is  to  be  displayed.  If  ISWITCH  is  one,  then  the  user  must  select  a  standard 
deviation  (SIGMA)  of  1.25  or  2.5  statute  miles.  An  ISWITCH  value  of  2  causes  an 
information  banner  to  be  displayed  and  the  value  of  SIGMA  is  set  to  0.5  statute  mile. 
ISWITCH  equal  to  3  indicates  that  the  user  has  chosen  to  enter  his  own  value  for  SIGMA. 
If  SIGMA  exceeds  2.5  statute  miles,  then  the  message  "Standard  Deviation  Out  of  Range" 
is  flashed  and  the  user  is  queried  to  enter  a  new  value.  All  values  of  SIGMA  are  converted 
to  feet  before  the  next  screen  is  displayed. 

Integer  variable  NUMAIR  counts  the  number  of  aircraft  configurations.  Its  veilue  is 
set  to  zero  prior  to  entering  loop  80.  Next,  NUMAIR  is  incremented  by  1  and  the  screen 
entitled  "Available  Aircraft  From  NOISEFILE"  is  displayed.  COMMON  block  NOISE 
passes  to  subroutine  OPS2  character  arrays  that  contain  key  descriptive  parameters  from 
NOISEFILE.  These  character  arrays  are  incrementally  displayed  when  the  arrow  or  page 
keys  are  pressed.  Briefly,  the  scrolling  feature  operates  by  having  the  keyboard 
continually  sampled  to  determine  if  a  key  has  been  pressed  and,  if  so,  whether  it  was  a 
key  that  constitutes  a  screen  action.  The  most  important  aspect  of  this  routine  is 
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recording  which  aircraft  file  in  the  NOISEFILE  data  base  has  been  selected.  Using  an 
integer  counter  variable  I,  the  NOISEFILE  recorded  number  is  stored  in  integer 
array  ICODE. 

When  drawing  the  operation  parameter  screen,  variable  ISWITCH  is  used  to  display 
the  heading  that  describes  the  mission  type.  Entry  labels  are  drawn  to  the  screen. 
If  integer  variable  INPUT  is  equal  to  1,  then  the  program  will  accept  a  centerline  track 
offset.  An  internal  read  and  a  print  window  library  call  are  used  to  draw  on  the  screen; 
the  aircraft  name,  the  power  setting,  and  the  air  speed.  Data  entry  is  achieved  by 
continually  sampling  the  keyboard,  checking  to  see  if  a  key  has  been  pressed  and,  if  so, 
whether  it  was  a  number,  a  RETURN,  or  an  ESC  key.  If  the  key  is  a  numerical  key, 
the  entry  is  read.  The  cursor  is  positioned  to  the  beginning  of  the  line;  the  entry  is 
then  echoed. 

Subroutine  OPS2  flashes  a  message  asking  if  another  configuration  is  desired.  If  the 
answer  is  yes,  then  program  execution  continues  at  loop  80.  If  the  answer  is  no, 
subroutine  OPS2  determines  whether  a  file  is  to  be  saved  on  disk  and  the  filename. 
Control  is  returned  to  the  main  program. 

5.3.7  OPS3 

This  routine  will  change  the  default  settings.  The  default  parameters  are  passed  to 
OPS3  through  COMMON  block  AIRPAR.  Entry  labels  are  drawn  and  the  current  default 
settings  are  written  to  the  screen.  The  keyboard  is  continually  sampled  to  determine  if  a 
key  has  been  pressed.  A  keystroke  will  cause  a  screen  action  only  if  an  ESC,  RETURN,  or 
numeric  key  is  pressed. 

5.3.8  CALC 

This  routine  calculates  the  noise  level  across  the  corridor.  There  are  two  parts  to 

this  program;  a  table  is  constructed  of  adjusted  mean  pressure  versus  offset  per  aircraft; 

in  the  second  part,  pressure  squared  weighted  by  the  probability  of  an  event  along  the 

track  per  aircraft  are  summed  for  each  ground  position.  Knowing  the  pressure  squared 

sum  across  the  track,  it  is  a  simple  calculation  to  determine  the  values  of  the 

L  .  ,  and  L  noise  metrics, 
dn’  eq 


32 


1 


The  table  construction  begins  with  an  array  of  25  sound  exposure  levels  (SEL) 
tabulated  as  a  function  of  profile  distance.  Using  these  SEL  values,  a  table  is  constructed 
for  the  slant  range  for  offset  values  ranging  from  0  to  73,000  feet  versus  SEL. 

The  pressure  squared  versus  offset  table  accounts  for  both  attenuations  and 
penalties.  Two  attenuation  mechanisms  considered  are  air  and  ground  absorption. 
OMEGA  10  uses  the  humidity  and  temperature  data  from  subroutine  OPS3  to  determine 
the  air  absorption  coefficient  and,  in  turn,  the  sound  exposure  level  reduction.  The  sound 
exposure  levels  stored  in  array  SEL  are  passed  to  CALC  via  COMMON  block  AIRPAR. 
Lateral  attenuation  is  accounted  in  subroutine  CALC  only  when  BETA,  the  elevation 
angle,  is  less  than  45  degrees.  This  equation  is  an  approximation  of  the  lateral 
attenuation  methodology  used  in  NOISEMAP  6.0  (Reference  6). 

Subroutine  CALC  accounts  for  increases  in  noise  intrusion  from  high-rate,  low- 
altitude  aircraft  by  calculating  the  Onset  rate  as  described  in  Reference  7.  The  exact 
form  of  this  equation  was  determined  by  curve  fitting  Onset  Rate  versus  Speed,  Altitude, 
Offset,  and  SEL  for  over  450  experimental  cases.  If  the  rate  is  greater  than  15  dB/sec 
and  less  than  30  dB/sec,  the  SEL  values  cire  penalized  by: 

Onset  Penalty  =  16.6  log^Q  (Onset/l5.0) 

If  the  onset  rate  is  greater  than  30  dB/sec,  then  a  5  dB  penalty  is  built  into  the 
tabulated  values. 

The  second  part  of  subroutine  CALC  determines  the  pressure  squared  along  the 
track,  accounting  for  span-wise  variability  in  aircraft  trajectory  and  the  different  types 
of  aircraft  flying  along  the  corridor.  The  algorithm  consists  of  three  nested  loops.  The 
innermost  loop  integrates  a  Gaussian  distribution  over  a  range  of  -2*SIGMA  -  10,000  feet 
to  2*SIGMA  +  10,000  feet,  while  the  second  loop  iterates  on  the  different  aircraft 
specified  by  the  subroutine  OPS2.  The  running  index  for  the  outer  loop  is  ground  position 
which  ranges  between  +50,000  feet. 

5.3.9  GAUSS(XINrr,XFINL,FMEAN,VAR,PROB) 

This  subroutine  computes  the  probability  over  an  interval.  XINIT  and  XFINL  are  the 
lower  and  upper  limits  in  the  integration.  VAR  is  the  variance  for  the  user-specified 
training  mission.  PROB  is  the  probability  computed  for  indicated  conditions. 
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5.3.10  FDMia,3,K,P) 

This  FUNCTION  uses  finite  difference  techniques  to  interpolate  the  tabulated  sound 
intensity  table  stored  in  FINTEN.  Array  FINTEN  is  passed  from  CALC  through  COMMON 
block  SPACE.  Four  integer  variables  are  passed  to  FINTEN  from  CALC;  variable  I 
specifies  the  configuration  number,  variable  J  points  to  the  present  corridor  position, 
variable  K  indicates  which  of  the  three  noise  metrics  is  to  be  used,  and  variable  P  is  used 
as  the  pivotal  point  for  the  interpolation. 

The  difference  equation  used  in  the  analysis  depends  on  pointer  J.  If  the  calculation 
is  at  the  top  or  bottom  of  the  difference  table,  the  equation  used  is  either  forward  or 
backward  technique,  respectively.  At  any  other  position  on  the  table,  the  technique 
employed  is  Stirling  central  difference. 

5.3.11  PLTMAIN(SWITCH,TOGGLE) 

PLTMAIN  is  the  main  plotting  program  for  ROUTEMAP.  Logical  variables  SWITCH 
and  TOGGLE  are  status  variables  passed  from  the  main  program.  SWITCH  equal  to  true 
indicates  that  the  plotting  data  resides  on  disk.  When  SWITCH  is  true,  PLTMAIN  queries 
the  user  for  a  filename.  A  "file  not  found"  message  is  displayed  only  once.  The  user  is 
given  two  tries  to  enter  the  correct  filename  before  PLTMAIN  returns  to  the  main 
program.  If  SWITCH  is  equal  to  false,  noise  level  data  is  passed  to  the  subroutine  through 
COMMON  block  RESULT. 

Subroutine  calls  are  made  to  PLTMENU  eind  DEVICE  to  determine  the  plot 
parameters  and  attributes.  These  parameters  in  turn  are  passed  between  plotting 
subroutines  through  COMMON  block  PLT.  The  FLAG  array  instructs  which  plots  the  user 
requests.  The  meaning  of  FLAG  equal  to  true  follows: 

•  f'LAG(l)  -  plot  Lj^^^ 

•  FLAG(2)  -  plot  Lj^ 

•  FLAG(3)  -  plot  L 

eq 

•  FLAG(4)  -  use  symbols  for  any  of  the  above  three  curves 

•  FLAG(5)  -  collapse  curves  on  to  one  plot 

•  FLAG(6)  -  plot  probability  of  high  annoyance 

Subroutine  PLTSPL  draws  the  noise  level  curves  if  any  one  of  the  first  three  FLAG 
elements  is  true.  A  true  value  for  FLAG(6)  enables  subroutine  PLTAPP  which  draws  the 
probability  of  high  annoyance  curve. 
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5.3.12  PLTMENU(SWrrCH,TOGGLE) 


This  subroutine  determines  the  plotting  parameters.  The  techniques  employed  in 
this  subroutine  are  analogous  to  previous  menus  that  require  data  entry.  Data  is  passed 
from  this  subroutine  through  the  PLT  COMMON  block.  Logical  variables  TOGGLE  and 
SWITCH  are  used  to  preserve  and  display  previously  entered  plot  menu  parameters  only 
when  the  sound  intensity  data  is  not  acquired  by  reading  a  file  from  the  disk. 

5.3.13  DEVICE 

This  subroutine  determines  the  lOPORT  and  MODEL  number.  A  screen  is  displayed 
offering  a  selection  from  five  possible  devices.  Variable  ICOUNT  records  the  option 
selection.  If  ICOUNT  is  less  than  or  equal  to  5,  then  data  arrays  lO  and  MOD  assign  a 
value  to  lOPORT  and  MODEL.  If  ICOUNT  equals  6,  then  a  screen  is  displayed  requesting 
the  user  to  input  the  model  number,  lOPORT  number,  and  scale  factor.  COMMON  block 
PLT  passes  the  parameters. 

5.3.14  PLTSPL 

This  subroutine  draws  the  noise  level  plots.  COMMON  blocks  RESULT  and  PLT  pass 
the  noise  levels  and  plotting  parameters  to  the  subroutine.  FMIN  and  FMAX  are  the 
corridor  coordinate  values  and  are  truncated  to  an  integer  multiple  of  500  feet.  NPTS  is 
the  number  of  points  used  to  plot  the  noise  level  curves. 

PLOT88  calls  are  used  extensively  throughout  the  subroutine.  These  subroutine  calls 
have  been  patterned  after  the  industry  standard  CALCOMP  style  calls  and  are  therefore 
interchangeable  with  other  similar  graphics  packages,  provided  adaptations  are  made  to 
ROUTEMAP's  graphics  calls.  The  program  initializes  PLOT88  using  the  lOPORT  and 
MODEL  number.  The  plot  is  scaled  to  correct  for  different  size  plotting  areas.  The 
program  centers  and  draws  the  plot  title  by  determining  the  number  of  characters  and, 
knowing  a  character  width,  measures  to  the  left  of  the  centerline  the  appropriate 
distance.  Logical  variable  FLAG  locations  I  through  3  determine  which  curves  are  to  be 
drawn.  FLAG  variables  4  and  5  are  used  to  determine  whether  symbols  are  to  be  drawn 
and  if  the  curves  are  to  be  collapsed  onto  one  plot. 
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5.3.15  PLTAPP 


This  subroutine  draws  the  probability  of  high  annoyance.  The  programming  logic 
used  for  plotting  parameter  initialization  and  the  drawing  of  axes  is  identical  to 
subroutine  PLTSPL.  DO  loop  10  contains  the  probability  calculations. 

5.3.16  DRAW(XARRAY,YARRAY,NPTS,FLAG1,INTEQ) 

This  subroutine  draws  a  line.  Variables  XARRAY  and  YARRAY  contain  the  absolute 
grid  coordinates  for  the  line.  NPTS  is  equal  to  the  number  of  line  segments  minus  one. 
A  true  value  for  FLAGl  instructs  the  subroutine  to  include  symbols  when  drawing  the  line. 
INTEQ  is  the  variable  defining  the  symbol  to  be  drawn. 

Stored  in  XARRAY  and  YARRAY,  position  NPTS  plus  one  is  the  minimum 
coordinate  value  and  in  position  NPTS  plus  two  is  the  separation  between  tic  marks. 
These  values  are  required  to  determine  the  X  and  Y  plotting  position.  The  program 
mechanics  is  based  on  a  series  of  move  and  draw  commands  using  PLOT 88  calls.  See 
PLOT88  reference  manual  for  detail  description  on  plotting  subroutine  library  calls. 


5.3.17  RANKODR 


This  subroutine  ranks  the  contributors  to  the  noise  exposure  at  a  specified  corridor 
position.  Logical  array  FLAG  is  assigned  a  true  value.  The  first  screen  requests  the 
distance  to  be  ranked  across  the  corridor.  A  message  is  flashed  if  the  distance,  variable 
DIST,  exceeds  the  corridor  limits  of  +50,000  feet.  Integer  variable  IV AL  points  to  the 


vector  in  the  RANK  matrix  whose  values  for  L  . 

dnmr 


correspond  to  the  distance  selected. 


COMMON  block  RESULT  passed  matrix  RANK  to  subroutine  RANKODR.  The  RANK 


matrix  contains  the  values  for  each  aircraft  type  at  500-foot  incremental  positions 

along  the  corridor. 


The  rank  contributors  screen  headings  are  drawn.  The  distance  value  selected  is 
echoed  on  the  screen.  Loop  80  contains  the  sorting  algorithm.  The  logic  to  this  algorithm 
is  a  derivative  of  a  bubble  sort,  the  difference  being  that  array  values  are  not  assigned; 
instead,  a  flag  is  switched  and  the  result  is  printed  on  the  screen.  A  FLAG  equal  to  true 


means  the  for  that  particular  aircraft  has  not  yet  been  identified  as  a  ranked 

element.  The  routine  begins  at  loop  50  where  it  looks  for  the  first  element  in  the  list  that 


has  not  been  ranked.  Once  it  is  found,  its  value  is  assigned  as  the  maximum  value, 
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VALMAX.  ITRACK  points  to  where  VALMAX  is  located  in  the  RANK  matrix.  Loop  60 
looks  through  the  RANK  vector  to  identify  a  noise  value  greater  than  VALMAX  that  has 
as  of  yet  not  been  listed  as  a  rank  element.  If  there  is  a  noisier  aircraft  that  meets  this 
condition,  then  both  VALMAX  and  ITRACK  are  reassigned.  Upon  the  completion  of  the 
DO  60  loop,  the  ITRACK  value  points  to  the  FLAG  which  is  set  to  false,  the  annoyance 
calculated,  and  the  results  printed  on  the  screen.  Loop  80  continues  to  iterate  until  all 
aircraft  are  ranked  and  in  turn  all  FLAGs  are  set  equal  to  false. 

5.3.18  ADDER 

This  subroutine  adds  the  noise  level  for  two  decibel  numbers.  Data  entry  is  achieved 
by  continually  sampling  the  keyboard,  checking  to  see  if  a  key  has  been  pressed,  and,  if  so, 
whether  the  key  was  a  number,  ESC,  or  RETURN  key.  If  the  key  is  a  numeric  key,  the 
entry  is  read  only  if  the  value  is  greater  than  or  equal  to  zero  or  less  than  or  equal  to 
1^0  dB.  When  both  the  entries  are  read,  the  sum  is  computed  and  echoed  beside  the 
equality  sign. 

5.3.19  FINISH 

This  routine  sets  the  computer  status  to  its  original  mode  prior  to  the  execution  of 
ROUTEMAP.  Displayed  on  the  terminal  is  the  screen  image  stored  in  INITSCN. 

5.3.20  BLOCK  DATA 

Relative  humidity,  air  temperature,  and  the  number  of  days  in  a  month  are  assigned 
in  the  data  block.  All  other  real  and  integer  COMMON  arrays  are  set  equad  to  a 
zero  value. 

5.4  COMMON  Blocks 

The  following  COMMON  blocks  and  their  contents  are  used  in  ROUTEMAP; 

•  /AIRPAR/  ALT(20),  FMEAN(20),  ICODE(20),  IREL,  ITEMP,  MONTH,  NUMAIR, 
FNUMD(20),  FNUMN(20),  POWER(20),  SEL(20,25),  SIGMA,  VEL(20) 

•  /FRAME/  SCN(1024) 

•  /NOISE/  MAX,  ACC(300),  AIRSPD(300),  DRAG(300),  NAME(300),  OPCR(300), 
POWDES(300),  POWSET(300),  POWUNT(300) 
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•  /PLT/  TITLE,  FMIN,  FMAX,  FLAG(6),  lOPORT,  MODEL,  FACT,  ICOUNT 

•  /RESULTS/  FLDNMR(201),  FLDNM(201),  FLEZ(201),  RANK(20,201) 

•  /SPACE/ FINTEN(20, 147,3) 


Figure  16  is  a  chart  showing  the  location  of  the  COMMON  blocks  in  the  programs. 
An  open  circle  indicates  the  COMMON  is  in  the  program.  A  solid  circle  indicates  that 
some  or  all  of  the  contents  of  the  block  are  defined  by  data  statements  within 
that  routine. 
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6.0  MAINTENANCE  MANUAL 

Program  maintenance  is  limited  to  modifying  ROUTEMAP's  input/output  (I/O) 
format  structure  to  reflect  any  format  I/O  change  in  OMEGA  10.  The  addition  of  new 
records  to  NOISEFILE  will  not  impact  ROUTEMAP  provided  the  format  structure  remains 
the  same  and  the  total  number  of  records  does  not  exceed  300.  Presently  NOISEFILE 
contains  243  records. 

6.1  OMEGA  10  Maintenance 

The  OMEGA  10  version  used  with  ROUTEMAP  has  been  modified  by  AAMRL  to 
accommodate  low-altitude  flights.  This  modification  resulted  in  an  output  file  that  had 
the  addition  of  two  more  cards  per  aircraft.  OMEGA  10  has  been  changed  to  facilitate 
data  entry.  The  following  items  document  those  changes: 

•  WRITE(1,9000),  WRITE(1,9100),  and  WRITE(1,9200)  are  now  comment  statements. 

•  READ(I,9300)SETUPFL,  READ{1,9300)  NOMAPFL,  and  READ(1, 9300)  NOISEFL 
are  also  comment  statements. 

•  T APE  5  filename  is  assigned  as  "OUT". 

•  T APE  3  filename  is  assigned  as  "IN". 

•  TAPE  7  filename  is  assigned  as  "NOISE". 

•  TAPE  6  filename  is  assigned  as  "OMEGA lO.OUT". 

This  completes  the  changes  made  to  OMEGA lO.FOR. 

6.2  ROUTEMAP  Maintenance 

If  the  I/O  of  OMEGA  10  were  modified,  then  ROUTEMAP  must  reflect  these 
changes.  The  input  deck  for  OMEGA  10  is  prepared  by  ROUTEMAP  in  the  MAIN  program 
beginning  at  line  154.  The  input  deck  is  written  to  TAPE  2  in  ROUTEMAP  having  OUT  as 
an  assigned  filename. 

The  output  deck  from  OMEGA  10  has  an  assigned  filename  of  IN  and  is  written  to 
disk.  Using  an  OPEN  statement  and  making  the  assignment  filename  equcU  to  IN  for 
TAPE  3,  ROUTEMAP  reads  OMEGA lO's  output  deck.  This  instruction  set  appears  in 
ROUTEMAP  beginning  at  line  199  in  the  MAIN  program. 
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6.3  Compiling  and  Linking 


To  compile  ROUTEMAP,  the  user  must  use  Microsoft  FORTRAN  Version  4.01  or 
later.  To  link  the  object  code,  PLOT88  and  SPINDRIFT  libraries  must  be  resident  on  the 
system.  The  command  procedure  used  to  compile  and  link  ROUTEMAP  follows: 

-  FL  /c  /Gt  MAIN.FOR 

-  FL  /c  /Gt  OPS.FOR 

-  FL  /c  /Gt  CALC.FOR 

-  FL  /c  /Gt  PLT.FOR 

-  FL  /c  /Gt  RANKODR.FOR 

-  FL  /c  /Gt  ADDER.FOR 

-  LINK/E  MAIN+OPS+CALC+PLT+RANKODR+ADDER  /SEG:1024,  ROUTEMAP,  , 
PLOT88  SPINLAB7 

This  procedure  is  embodied  in  the  batch  file  ROUTE.BAT. 

ROUTEMAP  can  be  run  on  a  system  with  only  512  kilobytes  of  memory  if  it  is  linked 
into  two  smaller  programs.  The  procedure  is  as  follows: 

-  LINK/E  MAIN+OPS+CALC+RANKODR+ADDER  /SEG:1024,  ROUTEl,  ,PLOT88 
SPINLAB7 

-  LINK/E  MAIN+PLT  /SEG:1024,  ROUTE2,  ,PLOT88  SPINLAB7 
T wo  error  messages  will  appear: 

•  Unresolved  externals: 

-  PLTMAI  in  file(s): 

MAIN.OBJ  (MAIN.FOR) 

•  Unresolved  externals: 

-  CALC  in  file(s): 

MAIN.OBJ  (MAIN.FOR) 

-  PAUSE  in  file(s): 

MAIN.OBJ  (MAIN.FOR) 

-  OPSl  in  file(s): 

MAIN.OBJ  (MAIN.FOR) 

-  ADDER  in  file(s): 

MAIN.OBJ  (MAIN.FOR) 
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-  OPS2  in  file(s): 

MAIN.OBJ  (MAIN.FOR) 

-  OPS3  in  file(s): 

MAIN.OB3  (MAIN.FOR) 

-  RANKODR  in  file(s): 

MAIN.OB3  (MAIN.FOR) 

To  execute,  the  user  types  ROUTE  1  then  proceeds  with  creating  a  configuration 
file,  calculating  the  noise  levels,  and  opening  a  file  to  store  the  results  on  disk.  Once  the 
file  has  been  created,  the  user  runs  ROUTE2  and  selects  the  plotting  option.  ROUTE2 
will  request  an  input  file  at  which  time  the  user  enters  the  filename  previously  specified. 

6A  Commercial  Software 

ROUTEMAP  requires  three  commercial  software  packages.  These  can  be  purchased 
from  the  manufacturers  or  from  dealers.  Care  must  be  taken  to  observe  license 
agreements  if  a  ROUTEMAP  load  module  is  transferred  to  a  different  computer. 

Microsoft  FORTRAN  is  manufactured  by  Microsoft,  16011  N.E.  36th  Way, 
Box  97017,  Redmond,  WA  98073-9717.  Software  support  is  provided  by  calling  206/ 
882-8089. 

The  SPINDRIFT  library  is  manufactured  by  Spindrift  Laboratories,  116  South  Har¬ 
vard  Avenue,  Arlington  Heights,  IL  60005-1644.  Software  support  is  provided  by  call¬ 
ing  312/255-6909. 

The  PLOT88  library  is  manufactured  by  PLOTWORKS,  Inc.,  P.O.  Box  12385, 
La  3olla,  CA  92037-0635.  Software  support  is  provided  by  calling  619/457-5090. 
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APPENDIX  A 


^  Computer  Program  Source  Listing 

C  ROUTEMAP 

C 

C  AUTHOR  :  MIKE  LUCAS 

C  SENIOR  ENGINEER 

C  WYLE  LABORATORIES 

C  2001  JEFFERSON  DAVIS  HWY. .  SUITE  701 

C  ARLINGTON.  VIRGINA  22202  (703)  892-6700 

C 

C  SPONSOR:  JERRY  D.  SPEARMAN 

C  BIODYNAMICS  ENVIRONMENT  BRANCH 

C  BIODYNAMICS  AND  BIOENGINEERING  DIV 

C  HARRY  G.  ARMSTRONG  AEROSPACE  MEDICAL 

C  RESEARCH  LABORATORY  (AFSC) 

C  WRIGHT-PATTERSON  AIR  FORCE  BASE.  OHIO  45433-6573 

C 

C  CONTRACT  NUMBER  F33615-85-3-0534  NOISEMAP  Forecast  Technology 

C 

c 

c 

c 

c 

c 

PROGRAM  MAIN 
SINCLUDE:  ‘ spinattr , hdr ‘ 

LOGICAL  CHECK.  FEXIST -  FLAG(6).  TOGGLE 
CHARACTER*!  IDUMMY 
CHARACTER»2  OPCR(300) 

CHARACTER*3  ACC(300),  AIRSPD(300) 

CHARACTER»5  POWSET(300) 

CHARACTER«6  POWUNT(300) 

CHARACTER*?  NAME(300) 

CHARACTER*! 6  FILNAM 
CHARACTER*20  POWDES(300) 

CHARACTER*25  DRAG (3 00) 

CHARACTER*80  INPUT.  OUTPUT.  TEXT(ll) 

JOB  CONTROL  CARD  FOR  OMEGAIO  INPUT.  ALL  BLANK  NO  PRINTED 
OUTPUT  AND  ONLY  SEL  VALUES  ARE  COMPUTED. 

COMMON  /AIRPARA/ALT(20)  .  FMEAN(20).  ICODE(20).  IREL .  ITEMP. 

1  MONTH.  NUMAIR,  FNUMD(20).  FNUMN(20).  POWER(20).  SEL(20,25). 

2  SIGMA,  VEL(20) 

COMMON  /FRAME/SCN( 1024 ) 

COMMON  /NOISE/MAX . ACC . AIRSPD . DRAG . NAME . OPCR . POWDES . POWSET . POWUNT 
COMMON  /RESULT/FLDNMR(201 ) .FLDNM(201 ) .FLEQ(201 ) .RANK(20 .201 ) 

DATA  TEXT(  1)/' RUNNING  OMEGAIO'/ 

DATA  TEXT(  2)/ 'POWER  OR  VELOCITY  OUT  OF  RANGE  FOR  OMEGAIO.  SEE  FIL 
IE  "OMEGAIO. OUT".  '/ 

DATA  TEXT(  3 )/' CALCULATING  NOISE  LEVELS'/ 

DATA  TEXT(  4)/ 'THE  INPUT  DATA'/ 

DATA  TEXT(  5)/'Config  Aircraft  DayOps  NiteOps  Altitude  Spe 
led  Power  Offset'/ 

DATA  TEXT(  6)/ 'ROUTEMAP  Version  1.0'/ 

DATA  TEXT(  7)/ 'TOTAL  Ldnmr  AT  DISTANCES  PERPENDICULAR  TO  ROUTE'/ 
DATA  TEXT(  8 ) / ' Dist ( Kf t )  dB  Dist(Kft)  dB  Dist(Kft)  dB  Dist(K 
1ft)  dB  Dist(Kft)  dB'/ 

DATA  TEXT  (  9)/ 'Save  Results  On  DisJc  (yes/no)  :'/ 

DATA  TEXT ( 10 ) / ' Fi le  already  exists.  Try  a  new  file  name.'/ 

A1 


c 

c 

c 

c 


n 

w 

C 


1000 

1100 

1200 

1300 


10 

20 


C 

C 


C 

C 

c 

c 

30 

C 

C 

40 


DATA  TEXT( 11 )/' Enter  File  Name  With  Extension  ;'/ 

INITIALIZE  WINDOW  ROUTINES 
CALL  INITIAL 

CHECK  THAT  NOISE  AND  OMEGA10.EXE  ARE  ON  THE  DISK 
IF(  .NOT.  FEXIST( 'NOISE- ) )  THEN 
CALL  FINISH 

WRITE(*.»)  'FILE  "NOISE"  NOT  FOUND' 

GO  TO  200 

ELSE 
END  IF 

IF(  .NOT.  FEXIST( 'OMEGA10.EXE' ) )  THEN 
CALL  FINISH 

WRITE(*.*)  'FILE  "OMEGAIO.EXE"  NOT  FOUND' 

GO  TO  200 
ELSE 
END  IF 
CALL  CLS 
CALL  CUROFF 

READ  NOISEFILE  DATA 

OPEN(UNIT=l .FILE*' NOISE' , STATUS* ‘ OLD ' .ERR=190) 

DO  10  I  =  1 .  300 

READd  .1000.END=20)  ACC(I).  OPCR(I) 

FORMAT (1 OX , A3 .A2 ) 

READd.  1100)  NAME(I).  AIRSPDd) 

FORMAT ( 3  7X . A7 . 1 OX . A3 ) 

READd.  1200)  DRAG(I) 

FORMAT ( 3 3X.A2S) 

READ (1.1300)  POWDES ( I )  . POWSET ( I ) . POWUNT ( I ) 
FORMAT (1 7X , A2  0 . IX . AS , IX . A6 ) 

READ  d . » ) 

READ (1 . « ) 

READ (1 . * ) 

CONTINUE 
CONTINUE 
CLOSE(UNIT»l ) 

MAX  =  I 

TRUE  TURNS  WINDOW  OPTION  —  ON  — 

FLAGd)  =  .TRUE. 

FLAG(2)  *  .TRUE. 

FLAG(3)  =  .FALSE. 

FLAG(4)  =  .TRUE. 

FLAG(5)  =  .TRUE. 

FLAG(6)  =  .FALSE. 

TOGGLE  *  .FALSE. 

DISPLAY  BANNER 
CALL  OPEN 

DRAW  MAIN  MENU 
CALL  MENU( IVAL .FLAG) 

GOTO  (40.  50.  60.  160.  170.  180.  190).  IVAL 

CONFIGURE  MILITARY  TRAINING  ROUTE 
IF(.NOT.  FLAGd))  GO  TO  30 
CALL  OPSl(ISWITCH) 

IFdSWITCH  .NE.  4)  CALL  OPS2(ISWITCH) 
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FLAG(l)  =  .FALSE. 

FLAG(2)  =  .TRUE. 

FLAG(3)  =  .TRUE. 

FLAG(4)  =  .TRUE. 

FLAG(5)  *  .TRUE. 

FLAG(6)  =  .FALSE. 

TOGGLE  =  .FALSE. 

GO  TO  30 
C 

C  SET  THE  TEMPERATURE.  HUMIDITY.  AND  NUMBER  OF  DAYS  IN  MONTH 
50  IF(.NOT.  FLAG(2))  GO  TO  30 

CALL  OPS3 
GO  TO  30 
C 

C  RUN  OMEGAIO  AND  CALC 

60  IF(.NOT.  FLAG(3))  GO  TO  30 

CALL  DEFWINDd  .8.0.14.79.SINGLBD.NORMAL) 

CALL  WINDBUF(SCN) 

CALL  CAPTSCN(SCN) 

CALL  CLRB(l) 

CALL  BLDBd  .2 .32  .TEXTd  )  .0.»»500) 

CALL  DISPSCN(SCN) 

CALL  WLOCATEd  .2 .32  .»500) 

C 

C  PREPARE  TEMPORARY  FILE  STRUCTURE 

IF(FEXIST( 'OUT'  )  )  CALL  ERASECOUT') 

IF(FEXIST(  ■  IN'  )  )  CALL  ERASECIN’) 

IF(FEXIST( 'OMEGAIO. OUT' ) )  CALL  ERASE (' OMEGAIO . OUT ' ) 

OPEN (UNIT=2 .FILE- 'OUT' .STATUS-'NEW' ,ERR-190) 

C 

C  PREPARE  INPUT  DECK  FOR  OMEGAIO 
C 

C  JOB  CONTROL  CARD  -  BLANK 

WRITE<2 .«) 

I  -  0 

70  I  -  I  +  1 

J  =  0 

80  J  =  J  +  1 

IF  (ACCdCODEd  )  )  .EQ.  ACC(  ICODE(  I+J )  )  .AND. 

1  I+J  .LE.  NUMAIR  .AND.  J  .LT.  12)  GO  TO  80 

C 

C  OUTPUT  PARAMETERS  FOR  EACH  AIRCRAFT  CARD  #1 

WRITE( 2 . 1400  )  ACC ( ICODE( I ) ) . ITEMP , IREL . J . POWUNT ( ICODE( I ) ) 
1400  FORMAT(A3 .13 .13 .10X.I2 .A6) 

C 

ASSIGN  1500  TO  IVAR 

IF  ( POWER  d)  .LT.  150)  ASSIGN  1600  TO  IVAR 
IF(POWER(I)  ,LT.  10)  ASSIGN  1700  TO  IVAR 
J  =  J  -  1 
IF(J  .GE.  6)  THEN 
WRITE(2.IVAR) 

1  (POWER(K) .IFIX(VEL(K) ) .OPCR( ICODE(K) ) .K-I+1 .K-I .1+5) 

WRITE(2 .IVAR) 

1  (POWER(K) .IFIX(VEL(K) ) ,OPCR( ICODE(K) ) . K-I +1 . K*I +6 . I +J ) 

ELSE 

WRITE(2.IVAR) 

1  (POWER(K) .IFIX(VEL(K) ) .OPCR( ICODE(K) ) .K-I+1 .K-I . I+J) 

END  IF 

1500  FORMAT(6(F5. 0.13 .A2 .12  )  ) 

1600  FORMAT(6(F5. 1 .13 .A2 .12) ) 
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1700 


FORMAT ( 6 ( F5 . 2 , 1 3 . A2 . 12 ) ) 

I  =  I  +  J 

IF(I  .LT.  NUMAIR)  GO  TO  70 
CLOSE(UNIT=2 ) 

C 

C  RUN  OMEGAIO 

CALL  EXEC (' OMEGAIO. EXE' . •  ■.«500) 

IF(  .NOT.  FEXIST( ‘ IN' ) )  THEN 
CALL  FINISH 

WRITE (*.*)  'FILE  "IN"  NOT  FOUND.' 

WRITE(*.*)  'OMEGAIO  FAILED  TO  WRITE  THE  FILE.' 

WRITE(*.«)  'CHECK  OMEGAIO. OUT  TO  DETERMINE  CAUSE  OF  ERROR.' 
GO  TO  200 
ELSE 
END  IF 
C 

C  READ  OMEGAIO  OUTPUT  DECK.  READ  INPUT  FORMATTED  FOR  OMEGAIO  VERSION 

C  THAT  BEGINS  CALCULATIONS  AT  AN  ALTITUDE  OF  100  FEET. 

OPEN(UNIT=3 .FILE=' IN' . STATUS= ' OLD ' .ERR=190) 

1  =  0 

90  READ( 3 .« ,END=100) 

1  =  1+1 

READ(3.1800,ERR=100)  (SEL(I.J).  J  =  1.  6) 

1800  FORMAT ( 22X .6F8. 1) 

READ ( 3 , * ) 

READ ( 3 , » ) 

READ ( 3 . »  ) 

READ(3 ,1900.ERR=100)  (SEL(I.J),  J  =  7.  14) 
READ(3.1900,ERR=100)  (SEL(I.J).  J  *  15,  22) 

1900  FORMAT ( 6X . 8F8 . 1 ) 

READ(3 ,2000,ERR*100)  (SEUI.J),  J  =  23.  25) 

2000  FORMAT(6X.3F8.1) 

READ ( 3 . «  ) 

READ ( 3 . * ) 

READ ( 3 . *  ? 

READ ( 3 . «  ) 

IF(I  .LT.  NUMAIR)  GO  TO  90 
CLOSE(UNIT=3 ) 

GO  TO  110 
C 

C  DISPLAY  SCREEN  :  POWER  OR  VELOCITY  OUT  OF  RANGE  FOR  OMEGAIO 

100  CLOSE(UNIT=3 ) 

CALL  WINDBUF(SCN) 

CALL  CLRB(l) 

CALL  BLDBCl ,2 .5.TEXT(2 ) .0,«500) 

CALL  DISPSCN(SCN) 

CALL  PAUSE 
FLAG(l)  =  .TRUE. 

FLAG(2)  =  .TRUE. 

FLAG(3)  =  .FALSE. 

FLAG(4)  =  .FALSE. 

FLAG(5)  =  .FALSE. 

FLAG(6)  =  .FALSE. 

TOGGLE  =  .FALSE. 

GO  TO  30 
C 

C  DISPLAY  SCREEN  :  RUNNING  CALC 

110  CALL  WINDBUF(SCN) 

CALL  CLRB(l) 

CALL  BLDBd  ,2.28.TEXT(3)  .  FLASH+NORMAL  .  *500  ) 
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CALL  DISPSCN(SCN) 

C 

C  EXECUTE  SUBROUTINE  CALC 

CALL  CALC 
C 

C  DISPLAY  INPUT  DATA 

CALL  DEFWIND(2 .0.0. 24.79, DBLBDR.NORMAL) 

CALL  WINDBUF(SCN) 

CALL  CLRB(2) 

I ADJUST  =  1 

IF(NUMAIR  .NE.  20)  THEN 

CALL  BLDB(2,  0.33, TEXT (  4),0.»500) 

I ADJUST  =  0 
ELSE 
END  IF 

CALL  BLDB(2.  3-IADJUST.  3.TEXT(  5).0.«500) 

CALL  DISPSCN(SCN) 

WRITE(OUTPUT,2100)  SIGMA/5280.0 
2100  FORMAT ( ‘ Standard  Deviation  *',F6.2,'  Statute  Miles') 

CALL  WLOCATE(2 .1-IADJUST.19.»500) 

CALL  PRINTW(2  .OUTPUTd  ;40)  ,»500) 

WRITEt OUTPUT. 2200)  ITEMP .  IREL.  MONTH 
2200  FORMAT ('Temp  =', 13 . 5X .' Humidity  »',I3.5X. 

1  ' Num  Days  In  Month  =  *,I2) 

CALL  WLOCATE (2,2-1 ADJUST . 1 3 . «500 ) 

CALL  PRINTW(2  .OUTPUTd  :54)  .*500) 

DO  120  1=1.  NUMAIR 

CALL  WLOCATE( 2 . I +3-1 ADJUST . 2 . »500 ) 

WRITE(OUTPUT.2300)  I  .NAME(  ICODE(  I )  )  .FNUMD(I )  .FNUMNd  )  . 

1  IFIX(ALT(I )  )  .IFIX(VEL(I) )  .POWER(I  )  .FMEANd  ) 

2300  FORMAT  (  3X,  12 .6X.A7.2X.F5. 1 ,4X  .F5. 1 , 5Xi'I5 , 6X  ,  13 .2X  .F6. 2  . 

1  3X.G8.3E1) 

CALL  PRINTW(2.0UTPUTd:73).«500) 

120  CONTINUE 

CALL  ANYKEY 
C 

C  DISPLAY  RESULTS  OF  CALCULATIONS 
CALL  WINDBUF(SCN) 

CALL  CLRB(2) 

CALL  BLDB(2.  0.30.TEXT(  6).0.*500) 

CALL  BLDB(2,  1,16. TEXT (  7),0,«500) 

CALL  BLDB(2.  2.  2. TEXT (  8).0.*500) 

CALL  DISPSCN(SCN) 

DO  130  I  =  1 ,  39 .  2 
IROW  =  (I+l )/2  +  2 
CALL  WLOCATE(2 ,IROW,2 .»500) 

WRITE( OUTPUT. 2400) 

1  ( ( J+I-101 )/2 ,IFIX(FLDNMR( J+I ) ) .  J  =  0.  160,  40) 

2400  FORMAT (2X. 13 .4X, 13 . 4 ( 5X . 13 . 4X . 13 ) ) 

CALL  FRINTW(2.0UTPUTd:74).«500) 

130  CONTINUE 

CALL  ANYKEY 
CALL  WINDBUF(SCN) 

CALL  CAPTSCN(SCN) 

CALL  CLRB( 1 ) 

CALL  BLDBCl .2.10, TEXT(  9).0,»500) 

CALL  DISPSCN(SCN) 

CALL  WLOCATEd  ,2 .41  ,*500) 

CALL  CURON 

READ (**,2500)  IDUMMY 
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2500  FORMAT (Al) 

IFdDUMMY  .EQ.  ‘V  .OR.  IDUMMY  .EQ.  'y')  THEN 
GO  TO  150 

140  CALL  WINDBUF(SCN) 

CALL  CLRB(l) 

CALL  BLDB( 1 . 2 . 10 .TEXT (10 ) .0.»500) 

CALL  DISPSCN(SCN) 

CALL  PAUSE 
CALL  CLS 

CALL  SYSTEM ( ‘DIR/W .*500) 

CALL  ANYKEY 

150  CALL  WINDBUF(SCN) 

CALL  CLRB(l) 

CALL  BLDBd  .2 .10  .TEXTdl )  .0.*500) 

CALL  DISPSCN(SCN) 

CALL  CURON 

CALL  WLOCATEd  .2 .42  .*500) 

READ(* .2600 .ERR=190  )  FILNAM 
2600  FORMAT (Al 6) 

CALL  CUROFF 

OPEN( UNIT=6 . ERR=140 . FILE=FILNAM . FORM= ' BINARY ' . RECL=1 ) 

WRITE(6)  FLDNMR.  FLDNM .  FLEQ 
CLOSE(UNIT=6) 

ELSE 
END  IF 

FLAG(l)  =  .TRUE. 

FLAG(2)  =  .FALSE. 

FLAG(3)  «  .FALSE. 

FLAG(4)  =  .TRUE. 

FLAG(5)  =  .FALSE. 

FLAG(6)  »  .TRUE. 

TOGGLE  =  .FALSE. 

GO  TO  30 
C 

C  PLOT  RESULTS 

160  CONTINUE 

C 

C  IF  OPS  IS  RUN  WITHOUT  RUNNING  CALC  THEN  DO  NOT  ALLOW  USER  TO  DRAW  PLOT 
C 

C  IF  FLAG(5)  TRUE  THEN  FILE  MUST  BE  ON  DISK. 

IF(.NOT.  FLAG(4))  GO  TO  30 

CALL  PLTMAIN(FLAG( 5) .TOGGLE) 

TOGGLE  =  .TRUE. 

GO  TO  30 

170  IF(.NOT.  FLAG(6))  GO  TO  30 

CALL  RANKODR 
GO  TO  30 

180  CALL  ADDER 

GO  TO  30 
C 

C  SPINDRIFT  ERROR 
500  CALL  FINISH 

IVAL  *  ISERROR( ) 

WRITE(*.*)  ‘SPINDRIFT  ERROR  NUMBER  -‘.IVAL 
GO  TO  200 
190  CALL  FINISH 

200  STOP  'Bye-Bye.  Have  a  Nice  Day  !’ 

END 

C 

C  A6 
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SUBROUTINE  INITIAL 

INTEGER«2  INITMOD,  VIDEOTYP.  IVIDMOD 
CHARACTER  VMODE«2 
LOGICAL  TESTMONO 
COMMON  /FRAME/SCN (1024) 

COMMON  /INTT/INITMOD.INITSCN 
DATA  TESTMONO/. FALSE. / 

DATA  IVIDMOD/3/ 

CALL  SPININIT 
INITMOD  =  VIDEOTYP() 

CALL  CAPTSCN(SCN) 

CALL  ATOS(SCN.INITSCN. 4096.0) 

CALL  NUMCLR 
CALL  INSCLR 

CALL  GETENVR( 'MONO' ,VM0DE.»10) 
TESTMONO  =  .TRUE. 

0  IF( INITMOD  .EQ.  7  .OR.  TESTMONO)  THEN 
RETURN 
ELSE 

CALL  VIDEOMOD( IVIDMOD) 

END  IF 
RETURN 
END 


SUBROUTINE  OPEN 

DRAWS  THE  OPENING  SCREEN  USING  A  ZOOM  FEATURE 


■  / 

Assessment  ‘ 


/ 


$INCLUDE:  ' spinattr . hdr ‘ 

CHARACTER»45  TEXT(8) 

COMMON  /FRAME/SCN (1024) 

DATA  TEXT(l)/'  ROUTEMAP 

DATA  TEXT(2)/'  Environmental  Noise 

DATA  TEXT(3)/'  ,  for'/ 

DATA  TEXT(4)/‘  Military  Aircraft  Training  Routes'/ 

DATA  TEXT(5)/'  by'/ 

DATA  TEXT(6)/'  United  States  Air  Force  AAMRL/BBE'/ 

DATA  TEXT( 7 ) / ' Wright-Patterson  Air  Force  Base,  Ohio  45433' 
DATA  TEXT(8)/'  Version  1.0'/ 


CALL  CUROFF 
DO  10  I  *  1.  12 

lULCR  =12-1 

lULCC  =  IFIX(39.48  -  FLOAT( I )«3 . 29 ) 

ILRCR  =12+1 

ILRCC  =  IFIX(39.52  +  FLOAT( I )»3 . 29 ) 

CALL  DEFWIND ( 1 . lULCR . lULCC . ILRCR . ILRCC . DBLBDR . NORMAL ) 
CALL  WINDBUF(SCN) 

CALL  CAPTSCN(SCN) 

CALL  CLRB(l) 

CALL  DISPSCN(SCN) 

CALL  TONE2( I«300+2000.I*300+1000,2 ) 
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10  CONTINUE 

CALL  WINDBUF(SCN) 

CALL  CLRB(l) 

CALL  BLDBd.  2 . 1 7  ,  TEXT ( 1  )  , BRIGHT  .  »20  ) 
CALL  BLDBd,  5 . 17  .TEXT(2  )  .BRIGHT  .«20  ) 
CALL  BLDBd.  7 . 17  .  TEXT  (  3  )  ,  BRIGHT  .  *20  ) 
CALL  BLDBd.  9 . 1 7  .  TEXT  (  4  )  .  BRIGHT  .  »2  0  ) 
CALL  BLDBd  ,  12 . 1 7  .  TEXT  (  5  )  .BRIGHT,  *20) 
CALL  BLDBd  .15.17.TEXT(6)  . BRIGHT. «20) 
CALL  BLDBd  ,16,17,  TEXT(  7)  .  BRIGHT.  »20) 
CALL  BLDBd  .  18 , 17  ,  TEXT  (  8  )  , BRIGHT,  »20) 
CALL  DISPSCN(SCN) 

CALL  WAITd  .27, *20) 

CALL  PAUSE 
20  RETURN 
END 


SUBROUTINE  MENU ( IVAL . FLAG ) 
INCLUDE:  ‘ spinattr . hdr ' 
INTEGER«2  KEYTYPE 
INTEGER  lAd) 

LOGICAL  ISKEY.FLAG(6) 
CHARACTER* 1  CHR,  GETCHX 
CHARACTER*60  TEXTdl) 
COMMON  /FRAME/SCNd024) 
DATA  TEXT (1)/* MAIN  MENU'/ 


DATA  TEXT(2)/' Please  Make  A  Selection'/ 


DATA  TEXT(3)/' (1 ) 
DATA  TEXT(4)/' (2) 
DATA  TEXT ( 5 ) /  '  ( 3 ) 
DATA  TEXT ( 6 ) / ' ( 4 ) 


Define  Operations  On  Military 
Change  Default  Settings.'/ 
Calculate  Noise  Level.'/ 

Plot  Noise  Level  Metrics.  Use 


Isk.  '  / 


Training  Route.'/ 


Data  Stored  On  Di 


DATA 
DATA 
1.  '/ 
DATA 
DATA 


TEXT ( 7 ) / ' ( 4 ) 
TEXT ( 8 ) / ' ( 5 ) 

TEXT ( 9 ) / ■ ( 6 ) 
TEXTdO)/'  (7) 


Plot  Noise  Level  Metrics.'/ 

Rank  Order  Contributors  To  Total  Noise  Exposure 

Decibel  Addition.'/ 

End  This  Program.'/ 


DATA  TEXTdl  )/■  Your  Selection 


0  CALL  DEFWINDd  .0,5,24, /-(j.DBLBDR. NORMAL) 
CALL  WINDBUF(SCN) 

CALL  CLS 

CALL  CAPTSCN(SCN) 

CALL  CLRBd) 

TITLE 

CALL  BLDBd.  2,31,  TEXT  (1  )  ,  0  .  *30 ) 

CALL  BLDBd,  3,24,  TEXT  ( 2  )  .  0  .  *30  ) 


OPERATIONS 
IF(FLAGd))  THEN 

CALL  BLDBd.  6.  8.  TEXT  (  3  ),  BRIGHT  .  *30  ) 
ELSE 

CALL  BLDBd.  6,  8,  TEXT  (  3  )  .  0  ,  *30  ) 

END  IF 


C 
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C  DEFAULT 

IF(FLAG(2))  THEN 


CALL 

ELSE 

BLDBd  . 

8. 

8. 

TEXT (4) .BRIGHT. *30) 

CALL 

BLDBd  . 

8. 

8. 

TEXT(4) .0.*30) 

END  IF 


CACULATE  NOISE  LEVEL 
IF(FLAG(3))  THEN 

CALL  BLDB( 1.10.  8.  TEXT ( 5 ) .BRIGHT . »30 ) 

ELSE 

CALL  BLDBd.lO.  8.  TEXT  (  5  )  .  0  .  »30  ) 

END  IF 

PLOT  MENU  ITEM 
IF(FLAG(4))  THEN 
IF(FLAG(5))  THEN 

CALL  BLDB(1.12.  8.  TEXT( 6 ). BRIGHT . »30 ) 

ELSE 

CALL  BLDB(1.12.  8.  TEXT( 7 ) .BRIGHT . »30 ) 

END  IF 
ELSE 

CALL  BLDB(1.12.  8.  TEXT ( 7 ) . 0 . *30 ) 

END  IF 

RANK  ORDER 
IF(FLAG(6))  THEN 

CALL  BLDB(1.14.  8.  TEXT ( 8 ) .BRIGHT , »30 ) 

ELSE 

CALL  BLDB(1.14.  8.  TEXT ( 8 ) . 0 . »30 ) 

END  IF 

END  PROGRAM 

CALL  BLDBd.ie.  8.  TEXT  (  9  ).  BRIGHT  .»30) 

CALL  BLDBd.18.  8.  TEXT  (10 )  .BRIGHT  .«30) 

CALL  BLDBd  .21 .40.TEXTdl )  .0.»30) 

CALL  DISPSCN(SCN) 

CALL  CURON 

CALL  WLOCATEd  .21 .56.W30) 

IF (.NOT.  ISKEY(CHR.KEYTYPE) )  THEN 
GO  TO  20 
ELSE 

CALL  CTOA(CHR.IA.O) 

IF((KEYTYPE  .EQ.  3)  .AND.  (IA(1).GE.49  .AND.  IA(1).LE.55))  THEN 
IVAL  =  lAd  )  -  48 
CHR  =  GETCH(KEYTYPE) 

ELSE 

CHR  =  GETCHX ( KEYTYPE ) 

CALL  WLOCATEd  .21 ,56.»30) 

CALL  PRINTWd  .  '  '  .*30) 

GO  TO  20 
END  IF 
END  IF 
CALL  CUROFF 
RETURN 
END 
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SUBROUTINE  FINISH 
INTEGER*2  INITMOD 
COMMON  /FRAME/SCN(1024) 
COMMON  /INIT/INITMOD.INITSCN 
CALL  VIDEOMOD (INITMOD) 

CALL  STOA(INITSCN.SCN,0,»10) 
CALL  DISPSCN(SCN) 

10  CALL  SPINWRAP 
RETURN 
END 


BLOCK  DATA 

COMMON  /AIRPARA/ALT(20)  .  FMEAN(20),  ICODE(20),  IREL,  ITEMP, 

1  MONTH.  NUMAIR,  FNUMD(20).  FNUMN(20).  POWER(20).  SEL(20.25). 

2  SIGMA,  VEL(20) 

COMMON  /RESULT/FLDNMR(201 ) ,FLDNM(201 ) .FLEQ(201 ) .RANK(20 .201  ) 

DATA  ALT/20*0. 0/.  FMEAN/20«0. 0/ .  ICODE/20»0/,  IREL/70/,  ITEMP/59/ 
DATA  MONTH/30/,  NUMAIR/O/,  FNUMD/20*0 . 0/ ,  FNUMN/20«0 . 0/ . RANK/4020* 
10.  0/ 

DATA  POWER/20»0. 0/ ,  SIGMA/0.0/,  VEL/20*0.0/ 

DATA  FLDNMR/201*0, 0/.  FLDNM/201«0 . 0/ .  FLEQ/201*0 . 0/ 

END 
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SUBROUTINE  OPSl ( ISWITCH  ) 

SINCLUDE:  ' spinattr . hdr ‘ 

INTEGER*2  KEYTYPE 
INTEGER  IA(1) 

LOGICAL  FLAG.  ISKEY 
CHARACTER*!  CHR .  GETCHX 
CHARACTER* 16  OLDNAM 
CHARACTER* 50  TEXT (10) 

COMMON  /AIRPARA/ALT(20) .  FMEAN(20).  ICODE(20).  IREL,  ITEMP. 

1  MONTH.  NUMAIR.  FNUMD(20).  FNUMN(20).  POWER(20).  SEL(20.25). 

2  SIGMA.  VEL(20) 

COMMON  /FRAME/SCN(1024) 

DATA  TEXTd)/' DEFINE  OPERATIONS  ON  MILITARY  TRAINING  ROUTE'/ 

DATA  TEXT(2 )/■ Select  A  Training  Mission'/ 

DATA  TEXT( 3 )/•(!)  Visual  Flying  Rules'/ 

DATA  TEXT(4)/'(2)  Instrument  Flying  Rules'/ 

DATA  TEXT(5)/' (3)  Special  Operations'/ 

DATA  TEXT(6)/‘(4)  Retrieve  Previously  Saved  Flight  Track  File'/ 
DATA  TEXT (7)/ 'Your  Selection  :'/ 

DATA  TEXT( 8 )/ ' Enter  File  Name  With  Extension  :'/ 

DATA  TEXT(9)/'File  Not  Found'/ 

DATA  TEXT (10)/ ‘Error  Reading  File'/ 

CALL  DEFWIND (1,0.5,15.73. DBLBDR . NORMAL ) 

10  CALL  WINDBUF(SCN) 

FLAG  *  .FALSE. 

CALL  CLS 

CALL  CAPTSCN(SCN) 

CALL  CLRB(l) 

CALL  BLDBd.  1 . 10  .  TEXT  (1 ) ,  0 ,  *500 ) 

CALL  BLDBd.  4 . 10  .  TEXT  ( 2  )  .  0  ,*500 ) 

CALL  BLDBd,  6 , 10  ,  TEXT  (  3  )  .  0  ,*500  ) 

CALL  BLDBd.  7 . 10  ,  TEXT(  4  )  ,  0  .*500 ) 

CALL  BLDBd,  8 , 1 0  .TEXT  (  5 )  ,  0  .*500 ) 

CALL  BLDBd.  9 . 10  ,  TEXT  (  6  ) ,  0  .*500 ) 

CALL  BLDB( 1,12.40, TEXT ( 7 ) . 0 , *500 ) 

CALL  DISPSCN(SCN) 

CALL  CURON 

20  CALL  WLOCATEd  .12 ,56,*500) 

IF (.NOT.  ISKEY(CHR.KEYTYPE) )  THEN 
GO  TO  20 
ELSE 

CALL  CTOA(CHR.IA.O) 

IF( (KEYTYPE  .EQ.  3)  .AND.  (IA(1).GE.49  .AND.  IA(1).LE.52))  THEN 
ISWITCH  =  lAd  )  -  48 
CHR  =  GETCHX (KEYTYPE) 

CALL  CUROFF 
ELSE 

CHR  *  GETCHX (KEYTYPE) 

CALL  WLOCATEd  .12 .56, *500) 

CALL  PRINTWd  .  '  '.*500) 

GO  TO  20 
END  IF 
END  IF 

CALL  DEFWIND (2.8,0.14.79, SINGLED . NORMAL ) 

IF( ISWITCH  .EO.  4)  THEN 


All 


GO  TO  40 

30  IF(FLAG)  GO  TO  10 

CALL  WINDBUF(SCN) 

CALL  CLRB(2  ) 

CALL  BLDB(2 .2 .10,TEXT(9) .0.«500) 

CALL  DISPSCN(SCN) 

CALL  PAUSE 
CALL  CLS 

CALL  SYSTEM ( -DIR/W  .»500) 

CALL  ANYKEY 
CALL  CLRKB 
FLAG  =  .TRUE. 

40  CALL  WINDBUF(SCN) 

CALL  CLRB(2) 

CALL  BLDB(2.2.10,TEXT(8) .0.*500) 

CALL  DISPSCN(SCN) 

CALL  CURON 

CALL  WLOCATE(2 .2 .42 .*500) 

READ(*.1100.ERR=40)  OLDNAM 
CALL  CUROFF 
1100  FORMAT (A1 6) 

OPEN( UNIT=4 . FILE=OLDNAM . STATUS* ‘ OLD ' . ERR*30 ) 
READ(4.1200.ERR=60)  IREL.  ITEMP.  MONTH.  NUMAIR.  SIGMA 
1200  F0RMAT(11X.I3. 3.9X.I3. 3 .13X.I2.2.17X.I2.2 .9X.F10. 1 ) 

DO  50  I  =  1 .  NUMAIR 

READ(4.1300.ERR-60)  ALT(I).  FMEAN(I).  FNUMD(I). 

1  FNUMN(I).  ICODE(I).  POWER(I).  VEL(I) 

1300  FORMAT(10X.F10.1.10X.F10.1.13X.F6.1.15X.F6.1,8X,I3.3.9X, 

1  F10.2.12X.F10.2) 

50  CONTINUE 

CL0SE(UNIT*4) 

ELSE 
ENDIF 
GO  TO  500 

60  CALL  WINDBUF(bCN) 

CALL  CLRB(2) 

CALL  BLDB(2.2.10.TEXT(10) .0**500) 

CALL  DISPSCN(SCN) 

CALL  PAUSE 
CL0SE(UNIT=4) 

GO  TO  10 
500  RETURN 
END 
C 

c 

c 

c 

c 

SUBROUTINE  OPS2 ( ISWITCH ) 

SINCLUDE:  ' spinattr . hdr ‘ 

DIMENSION  IA(1).  IROW(6).  ICOL(6).  VALUE(6) 

INTEGER*2  KEYTYPE 
LOGICAL  ISKEY 

CHARACTER  IDUMMY.  CHR.  GETCHX 
CHARACTER*2  OPCROOO) 

CHARACTER*3  ACC( 300 ) . AIRSPD( 300 ) 

CHARACTER»5  POWSET(300) 

CHARACTER*6  POWUNT(300) 

CHARACTER*?  NAME(300) 

CHARACTER* 16  NEWNAM 
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CHARACTER*20  POWDES(300) 

CHARACTER*25  DRAGOOO) 

CHARACTER«80  OUTPUT.  TEXT (50) 

COMMON  /AIRPARA/ALT(20)  .  FMEAN(20).  1CODE(20).  IREL,  ITEMP, 

1  MONTH.  NUMAIR.  FNUMD(20).  FNUMN(20).  POWER(20),  SEL(20.25). 

2  SIGMA.  VEL(20) 

COMMON  /FRAME/SCN(1024  ) 

COMMON  /NOISE/MAX , ACC . AIRSPD . DRAG . NAME . OPCR . POWDES . POWSET , POWUNT 
DATA  IROW  /8. 10. 12. 14. 16. 18/ 

DATA  ICOL  /32 .29.46,42 .45.47/ 

DATA  TEXTd)/' VISUAL  FLIGHT  RULES'/ 

DATA  TEXT  (2 )/ ‘  OPTION  (1)  -  Specific  well  definded  trac)c(s):  standa 
Ird'  / 

DATA  TEXT(3)/‘  deviation  is  1.25  statute  miles. '/ 

DATA  TEXT(4)/'  or'/ 

DATA  TEXT( 5 )/ ‘ OPTION  (2)  -  Multiple  scattered  tracks:  standard  dev 
liation ‘ / 

DATA  TEXT(6)/‘  is  2.5  statute  miles  about  centerline,' 

1/ 

DATA  TEXT(7)/'your  Selection  :'/ 

DATA  TEXT( 8)/' INSTRUMENT  FLYING  RULES'/ 

DATA  TEXT ( 9 )/' Assume  laterial  distribution  of  aircraft  relative'/ 
DATA  TEXT(IO)/'  to  route  centerline  described  by  a  Gaussian'/ 

DATA  TEXT(ll)/'  distribution  with  a  standard  deviation  of'/ 

DATA  TEXT(12)/'  0.5  Statute  Miles.'/ 

DATA  TEXT( 13 )/' SPECIAL  OPERATIONS'/ 

DATA  TEXT ( 14 )/' Enter  the  standard  deviation  in  Statute  Miles  (x.xx 
1)  :  '  / 


1 

1 


DATA 

DATA 

DATA 

DATA 


TEXT( 15 )/' Standard  Deviation  Out  of  Range'/ 
TEXT( 16 )/' Standard  Deviation  «'/ 

TEXT (17)/ 'AVAILABLE  AIRCRAFT  FROM  NOISEFILE'/ 


TEXT(18)/' 
DRAG 

DATA  TEXT(19)/' 
CONFIGURATION 
DATA  TEXT (20)/ 


AIR'/ 
AIRCRAFT 
SPEED ' / 


POWER 


DESCRIPTION 


POWER 


1 ' Use  arrow  keys  to  scroll  desired  aircraft  into  highlighted  area.' 
2/ 

DATA  TEXT(21 )/' press  RETURN  when  done.'/ 

DATA  TEXT (22 )/' OPERATION  PARAMETERS'/ 

DATA  TEXT (23)/' Power  (  )  ;'/ 

DATA  TEXT(24)/' Velocity  (  knots)  ;'/ 

DATA  TEXT(24)/' Speed  (  knots)  ;'/ 

DATA  TEXT( 25 ) / ‘ A1 titude  Above  Ground  Level  (feet)  :'/ 

DATA  TEXT( 26 )/■ Center 1 ine  Track  Offset  (feet)  :'/ 

DATA  TEXT(27 )/' Number  Of  Daytime  Monthly  Sorties  :'/ 

DATA  TEXT ( 28 )/' Number  Of  Nighttime  Monthly  Sorties  :'/ 

DATA  TEXT ( 29 )/■ press  ESC  when  all  done'/ 

DATA  TEXT ( 30 )/' Enter  Another  Configuration  (Yes/No)  :'/ 

DATA  TEXT(31 )/'Save  File  On  Disk  (Yes/No)  :'/ 

DATA  TEXT ( 32 )/' Enter  File  Name  With  Extension  :'/ 

DATA  TEXT( 33 )/ ' Fi le  already  exists.  Try  a  new  file  name.'/ 

GO  TO  (10.  30.  40).  ISWITCH 


VISUAL  FLIGHT  RULES  SCREEN 
10  CALL  DEFWINDd  .0.5.17.73.DBLBDR.NORMAL) 
CALL  WINDBUF(SCN) 

CALL  CLS 

CALL  CAPTSCN(SCN) 

CALL  CLRB(l) 
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CALL  BLDBd,  2 . 24  .  TEXT  ( 1  )  ,  0  .  »500  ) 

CALLBLDBd.  5.  7  .  TEXT  ( 2  )  .  0  .  *500  ) 

CALL  BLDBd.  6.  7  .  TEXT  (  3  )  ,  0  .  *500  ) 

CALLBLDBd,  8.  9  ,  TEXT  (  4  )  .  0  .  *500  ) 

CALL  BLDBd  ,10,  7  .  TEXT  (  5  )  ,  0  ,  *500  ) 

CALLBLDBd.il.  7  .  TEXT (  6  )  ,  0  ,  *500  ) 

CALL  BLDBd  .14,40,  TEXT(  7)  .0.*500) 

CALL  DISPSCN(SCN) 

CALL  CURON 

20  CALL  WLOCATEd  .14.56.*500) 

IF (.NOT.  ISKEY(CHR.KEYTYPE) )  THEN 
GO  TO  20 
ELSE 

CALL  CTOA(CHR.IA.O) 

IF((KEYTYPE  .EQ.  3)  .AND.  (IAd).GE.49  .AND.  IAd).LE.50))  THEN 
INPUT  =  lACl)  -  48 
CHR  =  GETCH ( KEYTYPE ) 

CALL  CUROFF 
ELSE 

CHR  =  GETCHX( KEYTYPE) 

CALL  WLOCATEd  .14. 56. *500) 

CALL  PRINTWd  .  '  ‘.*500) 

GO  TO  20 
END  IF 
END  IF 

IF( INPUT  ,EQ,  1)  THEN 
SIGMA  *  6600.0 

ELSE 

SIGMA  «  13200.0 

END  IF 
GO  TO  70 

INSTRUMENT  FLIGHT  RULES  SCREEN 
CALL  DEFWIND (2,0.5,11.73. DBLBDR , NORMAL ) 

CALL  WINDBUF(SCN) 

CALL  CLS 

CALL  CAPTSCN(SCN) 

CALL  CLRB(2) 

CALL  BLDB(2,  2,23.TEXT(  8).0.*500) 

CALL  BLDB(2.  5,10,TEXT(  9),0,*500) 

CALL  BLDB(2.  6 . 10 . TEXT ( 10 ) , 0 , *500 ) 

CALL  BLDB(2.  7 , 10 .  TEXTdl  ) .  0 .  *500 ) 

CALL  BLDB(2.  8 . 10 ,TEXT( 12 ) ,0 ,*500) 

CALL  DISPSCN(SCN) 

CALL  PAUSE 
SIGMA  =  2640.0 
INPUT  =  2 
GO  TO  70 

SPECIAL  OPERATIONS  SCREEN 
40  CONTINUE 

CALL  DEFWIND( 3 ,8. 0,14, 79, DBLBDR, NORMAL) 

CALL  WINDBUF(SCN) 

CALL  CLS 

CALL  CAPTSCN(SCN) 

CALL  CLRB(3) 

CALL  BLDB(3.2.5.TEXT(14) . NORMAL . *500 ) 

CALL  DISPSCN(SCN) 

CALL  CURON 
GO  TO  60 
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50  READ (*.1000)  IDUMMY 
1000  FORMAT (Al) 

CALL  WLOCATEO  .2 .58.*500) 

CALL  PRINTWO  .  ‘  *,*500) 

60  CALL  WLOCATEO  .2 ,58.*500) 

IF(.N0T.  ISKEy(CHR.KEYTYPE) )  THEN 
GO  TO  60 
ELSE 

CALL  CTOA(CHR.IA.O) 

IF((KEYTYPE  .EQ.  3)  .AND. 

1  ((IA(1).GE.46  .AND.  IA(1).LE.57)  .AND.  IA(1).NE.47))  THEN 

READ(* ,1100,ERR«50)  SIGMA 
1100  FORMAT(F6.2) 

ELSE 

CHR  =  GETCHX(KEYTYPE) 

CALL  WLOCATE( 3 .2 , 58.*500) 

CALL  PRINTW(3 . ■  *.*500) 

GO  TO  60 
END  IF 
END  IF 

IF(SIGMA  .GT.  2.5  .OR.  SIGMA  .LE.  0.0)  THEN 
CALL  CLRB(3) 

CALL  BLDB( 3,2.10, TEXT (15). NORMAL , *500 ) 

CALL  DISPSCN(SCN) 

CALL  CUROFF 
CALL  PAUSE 
GO  TO  40 
ELSE 
END  IF 
CALL  CLRB(3) 

CALL  BLDB (3.2.10. TEXT (16). NORMAL . *500 ) 

CALL  DISPSCN(SCN) 

CALL  WLOCATE(3 .2 .31 .*500) 

WRITE(».1200.ERR*30)  SIGMA 

1200  FORMAT ( F6 . 2 . '  Is  this  correct  (Yes/No)  :'.$) 

READ(*.1300)  IDUMMY 
1300  FORMAT (Al) 

IFdDUMMY  .EQ.  'N'  .OR.  IDUMMY  .EQ.  'n')  GO  TO  40 
SIGMA  =  SIGMA*5280.0 
INPUT  »  1 
CALL  CUROFF 
C 

C  BEGIN  AIRCRAFT  SELECTION 

70  NUMAIR  =  0 

CALL  DEFWIND (4.0.0.24,79, DBLBDR . NORMAL ) 

80  CONTINUE 

NUMAIR  =  NUMAIR  +  1 
I  =  -4 

90  CONTINUE 

CALL  WINDBUF(SCN) 

CALL  CLRB(4) 

CALL  BLDB(4,  1 . 23 . TEXT ( 17 ) . 0 . *500 ) 

CALL  BLDB(4.  3,  0 , TEXT ( 18 ) , 0 , *500 ) 

CALL  BLDB(4,  4.  0 , TEXT ( 19 ) , 0 . *500 ) 

CALL  BLDB(4.19.  8 . TEXT ( 20 ) . 0 , *500 ) 

CALL  BLDB(4.20.29.TEXT(21 ) ,0,*500) 

IF(I  .LT.-4)  I  »  MAX 

IF(I  .GT.  MAX)  I  -  -4 

C 

C  DISPLAY  TABLE  FROM  NOISEFILE  INCLUDE  SCROLLING  FEATURE 
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DO  100  J  =  7.  17 
lATTR  =  NORMAL 

IF(J  .EQ.  12)  lATTR  =  BRIGHT 

IF(I  .GE.  1  .AND.  I  .LE.  MAX)  THEN 

CALL  BLDB(4.J.  3 . NAME( I ) . lATTR . *500 ) 

CALL  BLDB( 4 . J . 12 . POWDES ( I ) . lATTR . »500 ) 

CALL  BLDB(4.  J.33  .POWSETd  )  .IATTR.»500) 

CALL  BLDB(4 , J.39.POWUNT(I ) .IATTR.»500) 

CALL  BLDB(4. J.47.DRAG(I ) ,IATTR.*500) 

CALL  BLDB( 4 . J , 73 . AIRSPD ( I ) . I ATTR . «500 ) 

ELSE 
END  IF 
1  =  1+1 
100  CONTINUE 

I  =  I  -  11 
CALL  DISPSCN(SCN) 

CHECK  TO  SEE  IF  THE  KEYBOARD  HAS  BEEN  PRESSED 
0  CONTINUE 

CHR  =  GETCHX(KEyTYPE) 

CALL  CTOA(CHR,IA.O) 

IF(KEYTYPE  .EQ.  4  .AND.  IA(1)  .EQ.  13)  GO  TO  120 
IF(KEYTYPE  .EQ.  8)  THEN 
IF(IA(1)  .EQ.  72)  THEN 
1  =  1-1 
GO  TO  90 
ELSE 
END  IF 

IF(IA(1)  .EQ.  73)  THEN 
I  »  I  -  10 
GO  TO  90 
ELSE 
END  IF 

IF(IA(1)  .EQ.  80)  THEN 
1  =  1+1 
GO  TO  90 
ELSE 
END  IF 

IF(IA(1)  .EQ.  81)  THEN 
I  =  I  +  10 
GO  TO  90 
ELSE 
END  IF 
ELSE 

GO  TO  110 
END  IF 

;0  CONTINUE 

RECORD  AIRCRAFT  CODE  FOR  OMEGAIO 
ICODE(NUMAIR)  =1+5 


SHOW  SCREEN  REQUESTING  DATA  INPUT 
CALL  DEFWIND(5. 0.5. 23 ,73 .DBLBDR. NORMAL) 
CALL  WINDBUF(SCN) 

CALL  CLS 

CALL  CAPTSCN(SCN) 

CALL  CLRS(5) 

IFCISWITCH  .EQ.  1)  THEN 

CALL  BLDB(5.  1 . 25 .TEXT ( 1 ) . 0 . »500 ) 
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c 

c 


130 

C 

C 


1400 

C 

c 


1500 

C 

C 


1600 

C 

C 


140 

1700 


150 


ELSE 

IFdSWITCH  .EQ.  2)  THEN 

CALL  BLDB(5.  1,23,  TEXT ( 8 ) . 0 . *500 ) 

ELSE 

CALL  BLDB(5.  1,26.  TEXT( 13 ) . 0 . »500 ) 

END  IF 
END  IF 

CALL  BLDB(5,  2 , 24 . TEXT ( 22 ) . 0 , »500 ) 

CALL  BLDB(5,  8 . 10 . TEXT ( 23 ) . 0 . «500 ) 

CALL  BLDB(5,10.10.TEXT(24)  .0,«500) 

CALL  BLDB(5.12.10.TEXT(25) ,0.»500) 

IFdNPUT  .EQ.  1)  CALL  BLDB(  5 . 14 , 10  .  TEXT  ( 26  )  .  0  ,  »500  ) 

CALL  BLDB(5,16.10,TEXT(27) .0.»500i 
CALL  BLDB(5.18,10,TEXT(28) .0.*500) 

CALL  BLDB(5.21 .22 .TEXT (29)  .0.»500) 

CALL  DISPSCN(SCN) 

CALL  CUEON 

RECCED  ZERO  VALUES  TO  ACCOUNT  FOR  WHEN  THERE  IS  NO  INPUT 
DO  130  1*1.5 
VALUEd)  =0.0 
CONTINUE 

WRITE  ARICRAFT  TYPE  AND  NUMBER  OF  CONFIGURATIONS  ON  SCREEN 
CALL  WLOCATE(5.5.10.»500) 

WRITE( OUTPUT.  1400)  NAMEdCODE( NUMAIR )  )  .  NUMAIR 

FORMAT ( ‘ Aircraft  Type  : ' , A7 , lOX , ' Conf iguration  :'.I2.$) 

CALL  PRINTW(5.0UTPUT(1 :49)  .»500) 

WRITE  POWER  ON  SCREEN 
CALL  WLOCATE(5.8,17.»500) 

WRITE( OUTPUT ( 1 : 12 ) . 1 500  )  POWSET ( ICODE ( NUMAIR ) ) . 

1  POWUNT dCODE( NUMAIR )  ) 

FORMAT (A5. IX. A6) 

CALL  PRINTW(5.0UTPUT(1 :12  )  .*500) 

WRITE  VELOCITY  ON  SCREEN 
CALL  WLOCATE(5.10.17.«500) 

WRITE(OUTPUT(l  :4)  .1600)  AIRSPDdCODE( NUMAIR )  ) 

CALL  PRINTW(5.0UTPUT(1 :4)  .«500) 

FORMAT (A3 .IX) 

INPUT  RESULTS.  ALLOW  FOR  EDITING  FEATURES 
1  =  1 
GO  TO  150 

READ (». 1700)  IDUMMY 
FORMAT ( A1  ) 

CALL  WLOCATE(5.IROW(  I  )  .ICOLd  )  ,*500) 

CALL  PRINTW(5.'  '.*500) 

CONTINUE 

IFdNPUT  .NE.  1  .AND.  I  .EQ.  4)  I  =5 
IFd  .EQ.  7)  I  =  1 

CALL  WLOCATE(5.IROW(I  )  ,  ICOLd  )  .»500) 

IF(.NOT.  ISKEY(CHR.KEYTYPE))  THEN 
GO  TO  150 
ELSE 

CALL  CTOA(CHR,IA.O) 

IF(KEYTYPE  .EQ.  16  .AND.  lA(l)  .EQ.  27)  GO  TO  170 
IF(KEYTYPE  .EQ.  4  .AND.  IA(1)  .EQ.  13)  THEN 
1  =  1+1 
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n  n  no 


1800 


160 

1900 


170 


2000 


CHR  =  GETCHX(KEYTYPE) 

GO  TO  150 
ELSE 
END  IF 

IF((KEYTYPE  .EQ.  3)  .AND. 

1  ((IA(1).GE.45  .AND.  IA(1).LE.57)  .AND.  IA(1).NE.47))  THEN 

CALL  PRINTW(5.‘  '.^SOO) 

CALL  WLOCATE ( 5 . IROW ( I) . ICOL ( I) , *500 ) 

READ(*,1800.ERR=140)  VAR 
FORMAT (GIO.  0) 

CALL  WLOCATEC  5 . IROW ( I ) . ICOL ( I) . *500 ) 

IF(I  .NE.  4  .AND.  VAR  .LT.  0.0)  GO  TO  160 
IF(I  .EQ.  4  .AND. 

1  (VAR  .LT.  -50000  .OR.  VAR  .GT.  50000))  GO  TO  160 

VALUE (I)  =  VAR 
WRITE(* . 1900 )  VALUE(I) 

FORMAT (GIO. 3 ,$ ) 

IF(VALUE(I)  .EQ.  VAR)  1=1+1 
GO  TO  150 
ELSE 

CHR  =  GETCHX(KEYTYPE) 

CALL  WLOCATE(5.IROW(I ) ,ICOL(I )  .*500) 

CALL  PRINTW(5. ■  ‘.*500) 

GO  TO  150 
END  IF 
END  IF 
CONTINUE 

CHR  =  GETCHX ( KEYTYPE ) 

COPY  INPUT  RESULTS  INTO  COMMON  ARRAYS 
POWER (NUMAIR)  =  VALUE(l) 

VEL(NUMAIR)  =  VALUE(2) 

ALT( NUMAIR)  *  VALUE(3) 

IF (INPUT  .EQ.  1)  THEN 

FMEAN( NUMAIR)  =  VALUE(4) 

ELSE 

FMEAN( NUMAIR)  =0.0 
END  IF 

FNUMD( NUMAIR)  =  VALUE(5) 

FNUMN( NUMAIR)  =  VALUE(6) 

CHECK  TO  SEE  IF  USER  WANTS  TO  ENTER  ANOTHER  AIRCRAFT 
CALL  DEFWIND( 6 .8.0,14, 79 . SINGLED . NORMAL ) 

CALL  CAPTSCN(SCN) 

CALL  WINDBUF(SCN) 

CALL  CLRB(6) 

CALL  BLDB(6,2 ,10,TEXT(30) ,0,*500) 

CALL  DISPSCN(SCN) 

CALL  CURON 

CALL  WLOCATE(6,2 ,48,*500) 

READ (*.2 000)  IDUMMY 
FORMAT ( A1 ) 

CALL  CUROFF 

IF(( IDUMMY  .EQ.  'Y'.OR.  IDUMMY  .EQ.  •y‘).AND.  NUMAIR  . LT .  20)  THEN 
CALL  CLS 
GO  TO  80 
ELSE 

CALL  WINDBUF(SCN) 

CALL  CLRB(6) 

CALL  BLDB(6,2 .10.TEXT(31 ) ,0,«500) 
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CALL  DISPSCN(SCN) 

CALL  CURON 

CALL  WLOCATE(6,2 ,38.*500) 

READ(«.2100)  IDUMMY 
2100  FORMAT (Al) 

IFdDUMMY  .EQ.  ‘Y'.OR.  IDUMMY  .EQ.  'y')  THEN 
GO  TO  190 

180  CALL  WINDBUF(SCN) 

CALL  CLRB(6) 

CALL  BLDB(6.2 ,10.TEXT(33) .0.»500) 

CALL  DISPSCN(SCN) 

CALL  PAUSE 
CALL  CLS 

CALL  SYSTEM ( 'DIR/W' .«500) 

CALL  ANYKEY 
CALL  CLRKB 

190  CALL  CAPTSCN(SCN) 

CALL  WINDBUF(SCN) 

CALL  CLRB(6) 

CALL  BLDB(6,2 .10.TEXT(32) .0.»500) 

CALL  DISPSCN(SCN) 

CALL  CURON 

CALL  WLOCATE(6,2.42.k500) 

READ(* .2200,ERR*180)  NEWNAM 
2200  FORMAT (Al 6) 

CALL  CUROFF 

OPEN ( UNIT*5 . FILE»NEWNAM . STATUS* ' NEW ' , ERR-1 80 ) 
WRITE(5.2300)  IREL,  ITEMP,  MONTH,  NUMAIR,  SIGMA 
2300  FORMATCRel.  Hum.  « *  .  13 . 3 .2X  ,  ‘  Temp.  » '  .  13 . 3  .2X  .  '  Num.  Days 
1  12. 2, 2X. ‘Num  Airplanes  * ’ , 12. 2 ,2X , ‘ Sigma  *‘,F10.1) 

DO  200  I  »  1,  NUMAIR 

WFITE(5,2400)  ALTd).  FMEAN(I).  FNUMD(I). 

1  FNUMN(I).  ICODE(I).  POWER(I).  VEL(I) 

2400  FORMAT ( ‘Altitude  * ‘ ,F10. 1 .2X , ‘ Of f set  «‘.F10.1,2X. 

1  Num.  Days  * ‘ , F6. 1 .2X , ‘ Num.  Nights  *‘.F6.1,2X, 

2  Code  =‘ ,13.3, 2X. ‘Power  = ‘ ,F10. 2 .2X , ‘ Vel ocity  =‘,F10.2) 

200  CONTINUE 

CLOSE(UNIT=5) 

ELSE 
END  IF 
END  IF 
500  RETURN 
END 
C 
C 
C 

c 

c 

SUBROUTINE  OPS3 
SINCLUDE:  ' spinattr . hdr ‘ 

DIMENSION  IA(1).  IROW(3),  ICOL(3).  IVALUEO) 

INTEGER*2  KEYTYPE 
LOGICAL  ISKEY 

CHARACTER  IDUMMY.  CHR ,  GETCHX 
CHARACTER*50  TEXT(5) 

COMMON  /AIRPAR/ALT(20)  .  FMEAN(20).  ICODE(20).  IREL,  ITEMP. 

1  MONTH.  NUMAIR.  FNUMD(20).  FNUM(20),  POWER(20).  SEL(20.25). 

2  SIGMA,  VEL (20) 

COMMON  /FRAME/SCN( 1024 ) 

DATA  TEXT(l)/‘ DEFAULT  SETTINGS ‘/ 
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n  o 


1000 

1100 

10 


20 

1200 


30 


1300 


DATA  TEXT ( 2 )/' Temperature  In  Degrees  Fahrenheit  :'/ 

DATA  TEXT( 3 )/' Relative  Humidity  : ' / 

DATA  TEXT( 4 )/‘ Number  Of  Days  In  Month  :’/ 

DATA  TEXT(5)/'press  ESC  when  all  done'/ 

DATA  IROW/3.6.9/ 

DATA  ICOL/45.29.35/ 

CALL  CURON 

CALL  DEFWIND (1,4,9,18,69. DBLBDR . NORMAL ) 

CALL  WINDBUF(SCN) 

CALL  CLS 

CALL  CAPTSCN(SCN) 

CALL  CLRB(l) 

CALL  BLDBd.  1  . 22  .  TEXT  ( 1)  .  0  .  *100  ) 

CALLBLDBCl,  3 . 10 , TEXT ( 2 ) , 0 . *100 ) 

CALL  BLDBd.  6 . 10  .  TEXT  (  3  )  .  0  .  *100  ) 

CALL  BLDBd  .  9 . 10  .  TEXT  ( 4  )  .  0  .  *100  ) 

CALL  BLDBd.l2.18.TEXT(5)  .0.*100) 

CALL  DISPSCN(SCN) 

IVALUEd)  =  ITEMP 
IVALUE(2)  =  IREL 
IVALUEO)  =  MONTH 
DO  10  I  »  1.  3 

CALL  WLOCATEd  .IROW(I)  .ICOL(I)  .*100) 

IF(I  .EQ.  1)  THEN 

WRITE(*.1000)  IVALUEd) 

ELSE 

WRITE(*.1100)  IVALUEd) 

END  IF 

FORMAT  dX,  13,  $) 

FORMATd3  .$  ) 

CONTINUE 

INPUT  RESULTS.  ALLOW  FOR  EDITING  FEATURES 
1  =  1 
GO  TO  30 

READ (*,1200)  IDUMMY 
FORMAT ( A1 ) 

CALL  WLOCATEd  .  IROWd  )  ,ICOL(I )  .*100) 

CALL  PRINTWd  .  '  '.*100) 

CONTINUE 

IFd  .EO.  4)  I  =  1 

CALL  WLOCATE( 1 ,IROW(I ) .ICOL(I ) ,*100) 

IF (.NOT.  ISKEY(CHR.KEYTYPE) )  THEN 
GO  TO  30 
ELSE 

CALL  CTOA(CHR,IA,0) 

IF(KEYTYPE  .EQ.  16  .AND.  IA(1)  .EQ,  27)  GO  TO  50 
IF(KEYTYPE  .EQ.  4  .AND.  IA(1 )  .EQ.  13)  THEN 
1  =  1+1 

CHR  =  GETCHX ( KEYTYPE ) 

GO  TO  30 
ELSE 
END  IF 

IF( (KEYTYPE  .EQ.  3)  .AND. 

1  ((IA(1).GE.45  .AND.  IA(1).LE.57)  .AND.  IA(1).NE.47))  THEN 

CALL  PRINTWd  ,  '  '  .*100) 

CALL  WLOCATEd  .IROWd  )  .ICOLd)  .*100) 

READ (*,1 300, ERR=20)  I VAR 
FORMAT (15) 

CALL  WLOCATEd  .IROWd  )  .ICOLd  )  .*100) 
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IF(I  .EQ.  1  .AND.  (IVAR.LT.-99  .OR.  IVAR. GT. 200 )  )  GO  TO  40 
IF(I  .EQ,  2  .AND.  (IVAR  .LT.  0  .OR.  IVAR  .GT.  100))  GO  TO  40 
IF(I  .EQ.  3  .AND.  (IVAR  .LT.  28  .OR.  IVAR  .GT.  31))  GO  TO  40 
IVALUE(I)  =  IVAR 

40  WRITE(»,1400)  IVALUE(I) 

1400  FORMAT ( 15. $) 

IFdVALUEd)  .EQ.  IVAR)  1  =  1+1 
GO  TO  30 
ELSE 

CHR  =  GETCHX ( KEYTYPE ) 

CALL  WLOCATEd  ,  IROW  ( I )  .  ICOL  (I )  ,«100) 

CALL  PRINTWd,'  '.»100) 

GO  TO  30 
END  IF 
END  IF 

50  CONTINUE 

CHR  =  GETCHX (KEYTYPE) 

CALL  CUROFF 
C 

C  COPY  INPUT  RESULTS  INTO  COMMON  ARRAYS 

ITEMP  =  IVALUE(l) 

IREL  =  IVALUE(2) 

MONTH  =  IVALUE(3) 

100  RETURN 
END 
C 
C 
C 
C 
C 

SUBROUTINE  PAUSE 
SINCLUDE:  ‘ spinattr . hdr ' 

COMMON  /FRAME/SCN(1024) 

CALL  DEFWIND(99,20.0.24.11 . DBLBDR . NORMAL ) 

CALL  WINDBUF(SCN) 

CALL  CLRB(99) 

CALL  BLDB(99.1.2.  'ANYKEY'  .NORMAL ,  **100  ) 

CALL  DISPSCN(SCN) 

CALL  ANYKEY 
100  RETURN 
END 
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SUBROUTINE  CALC 
LOGICAL  FEXIST 

COMMON  /AIRPARA/ALT(20) .  FMEAN(20),  ICODE(20).  IREL ,  ITEMP. 

1  MONTH.  NUMAIR.  FNUMD(20).  FNUMN(20).  POWER(20),  SEL(20,25). 

2  SIGMA.  VEL(20) 

COMMON  /RESULT/FLDNMR(201 ) .FLDNM(201 ) .FLE0(201 ) .RANK(20 .201 ) 
COMMON  /SPACE/FINTEN(20.147.3) 

DATA  PIDEG.DEGPI/0. 01745329252.57.29577951/ 

SLANT  RANGE  0  TO  73.000  FEET. 

DO  20  J  =  1 .  147 

OFFSET  =  FLOAT ( (J-1  )*500) 

DO  10  I  =1.  NUMAIR 

SLANT  =  S0RT( ALT(I)»ALT(I )  +  OFFSET » OFFS ET ) 

INT  =  IFIX( 10. 0*LOG10( SLANT)  -  19.0) 

IFdNT  .GT.  24)  INT  =  24 
IFCINT  .LT.  1)  INT  -  1 

51  =  10. 0««(FLOAT(INT+19)/10.  ) 

52  =  10. 0«»(FLOAT(INT+20)/10.  ) 

VALUE  =  SEL(I.INT)  + 

1  ( ( (SEL(I .INT)-SEL(I.INT+1 ) )/LOG10 ( S1/S2 )  ) 

2  *LOG10(SLANT/S1 ) ) 

BETA  =  ASIN(ALT(I)/SLANT)»DEGPI 

LATERAL  ATTENUATION 

IF(  BETA  .GT.  2.0  .AND.  BETA  .LT.  45.0)  THEN 
ATTEN  =  -0.3786  +  17.04/BETA 
ELSE 

IF(BETA  .LE.  2.0)  THEN 
ATTEN  *  8.14 
ELSE 

ATTEN  -  0.0 
END  IF 
END  IF 

VALUE  *  VALUE  -  ATTEN 
C 

C  EQUATION  FOR  PREDICTION  OF  ONSET  RATES  FOR  LOW-ALTITUDE. 

C  HIGH  SPEED  AIRCRAFT  OPERATIONS.  REFERENCE  LETTER  FROM 

C  JERRY  SPEAKMAN.  AAMRL  APAFB  FEBRUARY  9.  1987.  THIS 

C  EQUATION  APPEARS  IN  ZROUTE  FOR  THE  SURPRISE  EFFECT. 

C 

C  IF( OFFSET  .EQ.  0.)  OFFSET  -  1. 

C  ONSET  *  100.0/(1.0  +  2.7183*»(10.01  -  3 . 62*LOG10 ( VEL( I ) ) 

Cl  +  2.48»LOG10(ALT(I))  +  0 . 1 5«L0G1 0 ( OFFSET ) 

C  2  -  0. 0542«VALUE) ) 

C 

C  EQUATION  FOR  PREDICTION  OF  ONSET  RATE  RATE.  THIS  EQUATION 

C  WAS  REFERENCED  BY  JERRY  SPEAKMAN  ON  SEPTEMBER  18.  1987 

C 

ONSET  *  10. 0»«(0. 002247»VEL(I )  -  0 . 0003228*ALT ( I  ) 

1  -  0. 00014»OFFSET  t-  0 . 008987*VALUE  -  0.434) 

IF(ONSET  .LT.  15.)  THEN 
ONPEN  -  0.0 
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ELSE 

IF(ONSET  .GT.  30.)  THEN 
ONPEN  =  5.0 
ELSE 

ONPEN  =  16. 6«L0G10( ONSET/1 5.  0) 

END  IF 
END  IF 

INTENSITY  USED  FOR  Ldnmr 

FINTENd  .  J.l)  =  (  (FNUMD(I)+10.0»FNUMN(I)  )/ (FLOAT(  MONTH  )  *86400.  )  ) 
1  *10. »*( (VALUE  +  ONPEN)/10.0) 

INTENSITY  USED  FOR  Ldntn 

FINTENd  .J. 2  )  =  (  (FNUMDd  )+10.0*FNUMNd  )  )/( FLOAT ( MONTH  ) *86400 .  )  ) 
1  *10.**( (VALUE)/10.0) 

INTENSITY  USED  FOR  Leq 

FINTENd  .J.  3  )  =  (  (FNUMDd  )+FNUMNd  )  )/( FLOAT ( MONTH  ) *86400.  )  ) 

1  *10. **( (VALUE)/10. 0) 


0  CONTINUE 

0  CONTINUE 

TABLE  ALL  DONE  NOW  INTEGRATE  ACROSS  THE  MTR 
FIRST  DEFINE  SOME  CONSTANTS 

THE  CONSTANT  41  SERVES  AS  A  10.000  FOOT  BUFFER  ZONE.  THIS 
WILL  ENSURE  ALL  THE  NECESSARY  LEVELS  ARE  CALCULATED. 

KMAX  =  4*IFIX( SIGMA )/500  +  41 

IF (FLOAT ( KMAX ) /2. 0  .EQ.  FLOAT ( KMAX /2 ) )  KMAX  *  KMAX  +  1 
KARGU  »  (KMAX/2)  +  1 
DO  22  J  •  1.  201 
DO  23  I  =  1.  20 
RANKd.J)  »  0.0 
23  CONTINUE 

22  CONTINUE 
C 


c 

SIGMA 

KMAX 

KARGU 

c 

6600 

93 

47 

c 

13200 

147 

74 

c 

c 

2640 

63 

32 

c 

NOW  INTEGRATE 
DO  50  J  =  1. 

ACROSS 

201 

THE  ROUTE 

SUMl  =0.0 
SUM2  =0.0 
SUM3  =0.0 

GNDPOS  =  FLOAT ( (J-1 01 ) *500) 

DO  40  I  *  1 .  NUMAIR 

IF(GNDPOS.LT.FMEAN(I )-(2. 0*SIGMA+10000.  )  .OR. 

1  GNDPOS. GT. FMEANd  )  +  (2.0*SIGMA+10000.  )  )  GO  TO  40 

SUM4  =0.0 
DO  30  K  -  1.  KMAX 

DIFF  «  FLOAT ( (K-KARGU) *500. 0) 

X  =  FMEANd)  +  DIFF 
XINIT  =  X  -  249.0 
XFINL  «  X  +  250.0 

CALL  GAUSS (XINIT. XFINL. FMEAN (I ) . SIGMA. PROB) 
OFFSET  =  ABS( GNDPOS  -  X) 
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lOFF  =  IFIX(OFFSET/500. 0)  +  1 
P  =  MOD (OFFSET/500. 0,1.0) 

IF(P  .GT.  0.5)  THEN 
P  =  P  -  1 
lOFF  =  lOFF  +  1 
END  IF 

SUMl  »  SUMl  +  PROB*FDMI(I.IOFF.l.P) 

SUM2  =  SUM2  +  PR0B»FDMI(I,I0FF,2.P) 

SUM3  =  SUM3  +  PROB»FDMI(I .IOFF.3 ,P) 

SUM4  *  SUM4  +  PROB*FDMI(I .lOFF.l .P) 

30  CONTINUE 

IF(SUM4  .LE.  1.0)  THEN 
RANK(I.J)  =  0.0 
ELSE 

RANK(I.J)  =  10. 0*LOG10(SUM4 ) 

END  IF 

40  CONTINUE 

IF(SUM1  .LE.  1.0)  THEN 
FLDNMR(J)  =0.0 
ELSE 

FLDNMR(J)  =  10.0»LOG10(SUM1) 

END  IF 

IF(SUM2  .LE.  1.0)  THEN 
FLDNM(J)  *  0.0 
ELSE 

FLDNM(J)  =  10. 0»LOG10(SUM2) 

END  IF 

IF(SUM3  .LE.  1.0)  THEN 
FLEO(J)  =0.0 


ELSE 

FLEQ(J) 
END  IF 
CONTINUE 
RETURN 
END 


10. 0»LOG10(SUM3) 


FUNCTION  FDMI(I.J.K.P) 

FINITE  DIFFERENCE  METHOD  INTERPOLATION  FOR  NUMERICAL  FUNCTION 
TABULATED  AT  EQUAL  INTERVALS 


COMMON  /SPACE/FINTEN(20,147.3) 

IF(J  .EQ.  1)  THEN 
NEWTON-GREGORY  FORWARD  FORMULA 

FDMI  =  FINTENd  ,  J.K)  +  P» ( FINTEN(  I ,  J+1  .K )  -  FINTEN(  I  ,  J  , K )  ) 

1  +  (P«(P-1)/2.0)«(FINTEN(I.J+2.K)  -  2 . 0*FINTEN( I . J+1 . K ) 

2  +  FINTENd  .J.K)  ) 

ELSE 

IF(J  .EQ.  147)  THEN 
NEWTON-GREGORY  BACKWARD  FORMULA 

FDMI  =  FINTENd,  J.K)  +  P»(FINTENd  .  J  .K)  -  FINTENd  .J-1  .K)  ) 

1  +  (P»(P-l)/2.0)«(FINTENd.J.K)  -  2 . 0»FINTENd  .  J-1  .  K  ) 

2  +  FINTENd. J-2.K )  ) 

ELSE 

STIRLING  FORMULA 

FDMI  »  FINTENd.  J.K)  +  ( P/2 . 0  )«( FINTENd  .  J+1  .  K  )  -  FINTENd  ,  J-1  .  K  )  ) 
1  +  (P»P/2. 0)»(FINTENd  .J+1  .K)  -  2 . 0»FINTENd  .  J.K) 
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2  +  FINTENd  .J-1  .K)  ) 

END  IF 
END  IF 
RETURN 
END 


SUBROUTINE  GAUSS ( XINIT . XFINL . FMEAN . SIGMA . PROB ) 

DATA  PI/3.1415926536/ 

PARAMETERS  ARE  : 

XINIT  LOWER  LIMIT  FOR  X  VARIABLE 
XFINL  UPPER  LIMIT  ON  X  VARIABLE 
FMEAN  MEAN 

SIGMA  STANDARD  DEVIATION 

PROB  PROBABILITY  FOR  A  NORMAL  DISTRIBUTION 
Z1  *  (XINIT  -  FMEAN) /SIGMA 
Z2  =  (XFINL  -  FMEAN) /SIGMA 
CONST  »  1/(SQRT(2. 0»PI ) ) 

H  «  (Z2  -  Zl)/4.0 
SUM  =0.0 

IF((Z1  .LT.  -3.0  .AND.  Z2  .LT.  -3.0)  .OR.  (Z1  .GT.  3.0  .AND.  Z2 
1  .GT.  3.0))  GO  TO  20 
DO  10  J  *  1.  5 

X  »  FLOAT ( J-1 )»H  +  Zl 
P  *  EXP(-0. 5«X»«2  ) 

IF(J  .EQ.  1  .OR.  J  .EQ.  5)  THEN 
SUM  =  SUM  +  P 
ELSE 

IF(FLOAT(J)/2.0  .EQ.  FLOAT(J/2))  THEN 
SUM  «  SUM  +  4.0«P 
ELSE 

SUM  »  SUM  +  2.0«P 
END  IF 
END  IF 
10  CONTINUE 

20  PROB  =  (H/3. 0)*CONST«SUM 

RETURN 
END 
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MAIN  PLOTTING  PROGEAM 


SUBROUTINE  PLTMAIN ( SWITCH , TOGGLE ) 

INCLUDE:  ' spinattr . hdr ' 

INTEGER*2  lOPORT, MODEL 

LOGICAL  FLAG(6) .FLAGOl .SWITCH.TOGGLE 

CHARACTER  CHR 

CHARACTER»16  PTNAME 

CHARACTER»40  TITLE . TEXT ( 3 ) 

COMMON  /FRAME/SCN (1024) 

COMMON  /RESULT/FLDNMR(201 ) ,FLDNM(201 ) .FLEQ(201 ) ,RANK(20 ,201 ) 
COMMON  /PLT/TITLE . FMIN , FMAX . FLAG . lOPORT . MODEL . FACT . I COUNT 

DATA  TEXT(l)  /‘File  Not  Found'/ 

DATA  TEXT(2)  /‘Enter  Plot  File  Name  With  Extension  :'/ 

DATA  TEXT(3)  /'Error  Reading  File'/ 

IF (SWITCH)  THEN 
FLAGOl  *  .FALSE. 

IF  USER  DID  NOT  RUN  CALC  THEN  READ  IN  SPL  DATA 
CALL  DEFWINDd  .  8 , 0 . 14 , 79  .SINGLED , NORMAL  ) 

GO  TO  20 

10  IF(FLAGOl)  GO  TO  500 

CALL  WINDBUF(SCN) 

CALL  CLRB(l) 

CALL  BLDB(1.2.10,TEXT(1).0.»500) 

CALL  DISPSCN(SCN) 

CALL  PAUSE 
CALL  CLS 

CALL  SYSTEM( 'DIR/W' .»500) 

CALL  ANYKEY 
CALL  CLRKB 
FLAGOl  »  .TRUE. 

20  CALL  WINDBUF(SCN) 

CALL  CLS 

CALL  CAPTSCN(SCN) 

CALL  CLRB(l) 

CALL  BLDBd  ,2 .10.TEXT(2  )  .0.»500) 

CALL  DISPSCN(SCN) 

CALL  CURON 

CALL  WLOCATEd  ,2 .47,«500) 

READ (*.1000)  PTNAME 
1000  FORMAT (A1 6) 

CALL  CUROFF 

OPEN( UNIT=7 . ERR=10 . FILE=PTNAME. FORM- ' BINARY ' . RECL-1 . 

1  STATUS- 'OLD' ) 

READ( 7,ERR=30 )  FLDNMR.  FLDNM .  FLEQ 
CLOSE( UNIT-7) 

ELSE 
END  IF 
GO  TO  40 

30  CALL  WINDBUF(SCN) 

CALL  CLRBd) 

CALL  BLDBd  .2,10. TEXT(3)  ,0, *500) 

CALL  DISPSCN(SCN) 

CALL  PAUSE 


GO  TO  500 
40  CONTINUE 

C 

C  INPUT  PLOT 

CALL  PLTMENU( SWITCH .TOGGLE) 

C 

C  INPUT  THE  DEVICE 

CALL  DEVICE 
CALL  CLS 
C 

C  DRAW  THE  SOUND  PRESSURE  LEVELS 

IF(FLAG(1)  .OR.  (FLAG(2)  .OR.  FLAG(3)))  CALL  PLTSPL 
C 

C  DRAW  THE  PROBABILITY  OF  HIGHLY  ANNOYED  ACROSS  THE  MTR 

IF(FLAG(6))  CALL  PLTAPP 
C 

C  NOW  RETURN  TO  THE  MAIN  PROGRAM 

WRITE(».1100) 

1100  FORMAT (‘  READY  TO  RETURN  TO  MAIN  MENU.',/, 

1'  Please  press  <return>  to  continue. '.$ ) 

READ (». 1200)  CHR 
1200  FORMAT (Al) 

500  RETURN 

END 
C 
C 

C  DETERMINE  DEVICE  TYPE 

C 

c 

SUBROUTINE  DEVICE 
SINCLUDE:  ' spinattr . hdr ‘ 

DIMENSION  IA(1).  IROW(3),  ICOL(3) 

INTEGERW2  IO( 5 ) , lOPORT . IVALUE(2 )  .KEYTYPE .MOD( 5 ) .MODEL 

LOGICAL  FLAG ( 6 ) . ISKEY 

CHARACTER  IDUMMY . CHR . GETCHX 

CHARACTER»40  TITLE 

CHARACTER»50  TEXT (16) 

COMMON  /FRAME/SCN (1024) 

COMMON  /PLT/TITLE . FMIN . FMAX . FLAG , lOPORT . MODEL . FACT . ICOUNT 
C 

C  PLOTTING  DEVICE  DATA 

DATA  10  /  0.9600.  93.  97.  99/ 

DATA  MOD  /  5.  30.  93.  96.  99/ 

C 

C  POSITION  MARKERS  FOR  OTHER  PLOTTING  DEVICE  SCREEN 
DATA  IROW  /  8.  9.10/ 

DATA  ICOL  /35.34,34/ 

/-• 

DATA  TEXT(  1 )/' AVAILABLE  PLOTTING  DEVICES'/ 

DATA  TEXT(  2)/‘(l)  Epson  FX-80  Or  FX-85  Printer'/ 

DATA  TEXT(  3)/' (2)  HP  7475A  Plotter'/ 

DATA  TEXT (  4)/' (3)  Hercules  Graphics  Card  (HGC)'/ 

DATA  TEXT (  5)/' (4)  IBM  Enhanced  Graphics  Adapter  (EGA)'/ 

DATA  TEXT (  6)/' (5)  IBM  Color  Graphics  Adapter  (CGA)'/ 

DATA  TEXT (  7)/' (6)  Other'/ 

DATA  TEXT(  8)/' Your  Selection  :'/ 

DATA  TEXT(  9)/ 'OTHER  PLOTTING  DEVICE'/ 

DATA  TEXT( 10)/ 'Enter  the  lOPORT  number  and  MODEL  number.'/ 

DATA  TEXT ( 11 )/' Consul t  PLOT88  User  Reference  Manual  to  determine'/ 
DATA  TEXT ( 12 )/' appropriate  settings  for  your  device.'/ 
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DATA  TEXT(13 )/' lOPOKT  number  : ' / 

DATA  TEXT( 14)/ 'MODEL  number  :'/ 

DATA  TEXT( 15)/' Scale  Factor  ;'/ 

DATA  TEXT ( 16 )/' press  ESC  when  all  done'/ 

DETERMINE  THE  DEVICE 

CALL  DEFWINDd  .0,5,20.73  .DBLBDR, NORMAL) 

CALL  WINDBUF(SCN) 

CALL  CLS 

CALL  CAPTSCN(SCN) 

CALL  CLRB(l) 

CALLBLDBd,  1  . 20  .  TEXT  (1  )  ,  0 ,  »500 ) 

CALL  BLDBd,  4 . 12  .  TEXT  ( 2  )  .  0 .  «500  ) 

CALLBLDBd.  6 . 12  .  TEXT  (  3  )  .  0  .  «500 ) 

CALLBLDBd.  8 . 12  .  TEXT  ( 4  )  .  0 .  *500  ) 

CALL  BLDBd  .  10 , 12  .  TEXT  (  5  )  ,0.«500) 

CALL  BLDBd  .  12 . 12  .  TEXT  (  6  )  ,0.»500) 

CALL  BLDB( 1,14.12, TEXT ( 7 ) , 0 . *500 ) 

CALL  BLDBd  .17.40,TEXT(8)  .0.*500) 

CALL  DISPSCN(SCN) 

CALL  CURON 

0  CALL  WLOCATEd  ,17,56.*500) 

IF (.NOT.  ISKEY(CHR.KEYTYPE) )  THEN 
GO  TO  10 
ELSE 

CALL  CTOA(CHR.IA.O) 

IF((KEYTYPE  .EQ.  3)  .AND.  (IAd).GE.49  .AND.  IAd).LE.54))  THEN 
ICOUNT  »  lAd)  -  48 
CHR  «  GETCHX(KEYTYPE) 

ELSE 

CALL  WLOCATEd  ,17. 56. *500) 

CALL  PRINTWd  .  ‘  *,*500) 

CHR  *  GETCHX(KEYTYPE) 

GO  TO  10 
END  IF 
END  IF 

IFdCOUNT  .LE.  5)  THEN 
CALL  CUROFF 
lOPORT  «  10 (ICOUNT) 

MODEL  *  MOD (ICOUNT) 

ELSE 

CALL  DEFWIND (2.0,5,15.73, DBLBDR . NORMAL ) 

CALL  WINDBUF(SCN) 

CALL  CLS 

CALL  CAPTSCN(SCN) 

CALL  CLRB(2) 

CALL  BLDB(2.  1,23.TEXT(  9),0.*500) 

CALL  BLDB(2.  4 . 13  ,  TEXTdO  )  ,  0  ,*500  ) 

CALL  BLDB(2,  6,  9 . TEXT ( 11 ) . 0 . *500 ) 

CALL  BLDB(2,  7 . 15 . TEXT ( 12 ) . 0 . *500 ) 

CALL  BLDB(2.  8 . 20 .TEXT(13 ) . 0 , *500 ) 

CALL  BLDB(2,  9 , 20 . TEXT (14 ) , 0 . *500 ) 

CALL  BLDB(2,10,20,TEXT(15) .0.*500) 

CALL  BLDB(2 .13 .22 , TEXT (16 ) . 0 . *500 ) 

CALL  DISPSCN(SCN) 

INPUT  RESULTS  ALLOW  FOR  EDITING  FEATURES 
I  »  1 
GO  TO  30 

20  READ (*,1000)  IDUMMY 
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CALL  WLOCATE<2 , IROW ( I ) , ICOL ( I ) . *500 ) 

CALL  PRINTW(2 . •  -.*500) 

30  CONTINUE 

IF(I  .EQ.  4)  I  =1 

CALL  WLOCATE(2  .  IROW(I)  .ICOLd  )  .*500) 

IF (.NOT.  ISKEY(CHR.KEYTYPE) )  THEN 
GO  TO  30 
ELSE 

CALL  CTOA(CHR.IA.O) 

IF(KEYTYPE  .EQ.  16  .AND.  IA(1)  .EQ.  27)  GO  TO  40 
IFCKEYTYPE  .EQ.  4  .AND.  IA(1)  .EQ.  13)  THEN 
1  =  1+1 

CHR  =  GETCHX(KEYTYPE) 

GO  TO  30 
ELSE 
END  IF 

IF((KEYTYPE  .EQ.  3)  .AND. 

1  ( (IA(1 ).GE.46  .AND.  IA(1).LE.57)  .AND.  IA( 1 ) . NE. 47 ) )  THEN 

IF(I  .NE.  3)  THEN 

CALL  PRINTW(2 . ‘  -.wSOO) 

CALL  WLOCATE(2.IROW(I).ICOL(I) .*500) 
READ(*.1100.ERR«20)  IVALUE(I) 

CALL  WLOCATE(2.IROW(I).ICOL(I) .*500) 

WRITE(*.1200)  IVALUE(I) 

1100  FORMAT (15) 

1200  FORMAT(I5,$) 

1  =  1+1 
GO  TO  30 
ELSE 

CALL  PRINTW(2.*  ‘.*500) 

CALL  WLOCATE(  2  .  IRC«4  ( I )  .  ICOL  ( I  )  ,  *500  ) 

READ(*.1300,ERR»20)  FACT 

CALL  WLOCATE(2  .IR<0W(I )  .ICOLd  )  .*500) 

WRITE(*.1400)  FACT 
1300  FORMAT(F6.2) 

1400  FORMAT ( F6.2 .$ ) 

1  =  1+1 
GO  TO  30 
END  IF 
ELSE 

CHR  =  GETCHX(KEYTYPE) 

CALL  WLOCATE(2  .IROWd  )  .ICOLd  )  .*500) 

CALL  PRINTW(2 . '  '.*500) 

GO  TO  30 
END  IF 
END  IF 

40  CONTINUE 

CHR  =  GETCHX ( KEYTYPE ) 

CALL  CUROFF 
lOPORT  =  I’ALUEd) 

MODEL  =  IVALUE(2) 

END  IF 

CALL  WLOCATEd  .0,0.*500) 

500  RETURN 
END 
C 

c 

C  LINE  DRAWING  ROUTINE 

C 
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SUBROUTINE  DRAW ( XARRAY . YARRAY . NPTS . FLAGl . INTEQ ) 

DIMENSION  XARRAY(203).  YARRAY (203) 

INTEGER*2  NPTS 
LOGICAL  FLAGl .  FLAG2 
CHARACTER  CTEXT 
XMIN  =  XARRAY ( NPTS+1 ) 

YMIN  =  YARRAY ( NPTS+1 ) 

XDELTA  =  XARRAY (NPTS +2 ) 

YDELTA  =  YARRAY (NPTS +2) 

FLAG2  *  .TRUE. 

IF (FLAGl)  THEN 

IFdNTEO  .EO.  1)  IVAL  =  11 
IF( INTEQ  .EQ.  2)  IVAL  =  9 
IF( INTEQ  .EQ.  5)  IVAL  =  7 
ELSE 
END  IF 

DO  10  I  =  1,  NPTS 

XPOS  =  (XARRAY(I)  -  XMIN) /XDELTA 
YPOS  =  (YARRAY(I)  -  YMIN) /YDELTA 
IF(YPOS  .LE.  0.0)  THEN 
IF(FLAG2)  THEN 

CALL  PLOT(XPOS.O. 0.3) 

ELSE 

CALL  PLOT(XPOS.0.0.2) 

END  IF 

FLAG2  -  .TRUE. 

ELSE 

FLAG2  »  .FALSE. 

IF (FLAGl  .AND.  FLOAT ( I )/FLOAT( IVAL ). EQ. FLOAT ( I/I VAL ) )  THEN 
CTEXT  «  CHAR (INTEQ) 

CALL  SYMBOL ( XPOS , YPOS ,0.12. CTEXT . 0 . 0 . -2 ) 

ELSE 

IF(I  .EQ.  1)  THEN 

CALL  PLOT (XPOS. YPOS. 3) 

ELSE 

CALL  PLOT (XPOS. YPOS. 2) 

END  IF 
END  IF 
END  IF 

10  CONTINUE 
RETURN 
END 


PLOTTING  MENU 


SUBROUT I NE  PLTMENU ( SW I TCH . TOGGLE ) 

SINCLUDE:  ‘ spinattr . hdr ’ 

DIMENSION  IA(1).  IROW(9),  ICOL(9).  VALUE(2) 

INTEGER«2  lOPORT . KEYTYPE , MODEL 
LOGICAL  FLAG(6) .ISKEY, SWITCH. TOGGLE 
CHARACTER  CHR .  GETCHX .  CTEXT 
CHARACTEE»3  IDUMMY 
CHARACTER»40  TITLE 
CHARACTER*65  TEXT (17) 

COMMON  /FRAME/SCN (1024) 

COMMON  /PLT /TITLE . FMIN . FMAX . FLAG . lOPORT , MODEL . FACT . I COUNT 
^  A30 
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DATA  IROW/  4.  7.  8,11.12.13.14,16.20/ 

DATA  ICOL/23 .29,29,46,44 .44 ,45,36.37/ 

DATA  TEXT(  l)/'PLOT  RESULTS'/ 

DATA  TEXT (  2 ) / ' GENERAL ' / 

DATA  TEXT(  3)/ 'Enter  Title  :'/ 

DATA  TEXT(  4)/'Enter  The  Distances  Across  MTR  To  Be  Plotted.'/ 
DATA  TEXT (  5)/ 'Minimum  :'/ 

DATA  TEXT (  6)/ 'Maximum  :'/ 

DATA  TEXT(  7)/ 'NOISE  LEVEL'/ 

DATA  TEXT(  8)/ 'Do  You  Want  To  Plot  Ldnmr  (Yes/No)  :'/ 

DATA  TEXT(  9)/ 'Do  You  Want  To  Plot  Ldn  (Yes/No)  :'/ 

DATA  TEXT(10)/'Do  You  Want  To  Plot  Leq  (Yes/No)  :'/ 

DATA  TEXT(ll)/'Do  You  Want  Symbols  Used  (Yes/No)  :'/ 

DATA  TEXT(12)/'Do  You  Want  To  Collapse  These  Curves'/ 

DATA  TEXT(13)/'Onto  One  Plot  (Yes/No)  :'/ 

DATA  TEXT (14)/ 'ANNOYANCE'/ 

DATA  TEXT(15)/'Do  You  Want  To  Plot  The  Probability  Of'/ 

DATA  TEXT(16)/'High  Annoyance  (Yes/No)  :'/ 

DATA  TEXT( 17)/ 'press  ESC  when  all  done'/ 

DRAW  PLOTTING  MENU 

CALL  DEFWINDd  .0.5,24.73  .  DBLBDR .  NORMAL  ) 

CALL  WINDBUF(SCN) 

CALL  CLS 

CALL  CAPTSCN(SCN) 

CALL  CLRBd  ) 

CALL  BLDBd,  1,28.TEXT(  1),0,«500) 

CALL  BLDBd.  3,  5. TEXT (  2).0.»500) 

CALL  BLDBd.  4.10.TEXT(  3).0.»500) 

CALL  BLDBd,  6.10.TEXT(  4).0.»500) 

CALL  BLDBd.  7. 20. TEXT (  5).0.»500) 

CALL  BLDBd.  8,20.TEXT(  6).0,«500) 

CALL  BLDBd, 10.  5.TEXT(  7).0,»500) 

CALL  BLDBd  .11 ,10.TEXT(  8).0.»500) 

CALL  BLDBd  ,12 .10.TEXT(  9).0.«500) 

CALL  BLDBd  .13,10.TEXTd0)  ,0.«500) 

CALL  BLDB (1,14,10. TEXT ( 1 1 ) . 0 . »500 ) 

CALL  BLDBd  .15.10.TEXTd2  )  ,0.«500) 

CALL  BLDBd  .16.12.TEXT(13)  .0.K500) 

CALL  BLDBd, 18,  5  .  TEXT  ( 14  ) ,  0  ,  *500 ) 

CALL  BLDBd  ,19.10.TEXT(15)  ,0.»500) 

CALL  BLDBd  .20.12  ,TEXT(  16  )  ,  0  .«500  ) 

CALL  BLDBd  ,  22 . 22  .  TEXT  ( 17  )  ,0.»500) 

CALL  DISPSCN(SCN) 

CALL  CURON 

IF (SWITCH  .OR.  .NOT.  TOGGLE)  THEN 

WRITE  DEFAULT  VALUES  TO  SCREEN 

CALL  WLOCATEd  .  IROW(2  ) ,  ICOL(2  ) ,  *500  ) 

CALL  PRINTWd  .' -50000.  '  .«500) 

CALL  WLOCATEd  .IROW (3  )  ,ICOL(3)  .«500) 

CALL  PRINTWd  ,' 50000.  '  ,»500) 

INPUT  RESULTS,  ALLOW  FOR  EDITING  FEATURES 
VALUE( 1 )  =  -50000. 0 
VALUE(2)  =  50000.0 

ELSE 

CALL  WLOCATE( 1 , IROW(l ) .ICOL( 1 ) .»500) 

WRITE(».1200)  TITLE 

CALL  WLOCATEd  .IROW(2)  .ICOL(2)  ,»500) 
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WRITE(« . 1400)  VALUE(l) 

CALL  WLOCATEd  .IR0W(3  )  .ICOLO  )  .«500) 

WRITE(* .1400)  VALUE(2) 

DO  10  I  =  1  ,  6 

CALL  WLOCATEd  ,  IROW(I+3  )  .ICOLd+3  )  .*500  ) 

IF(FLAGd))  THEN 

CALL  PRINTWd  .•  Yes'  .*500) 

ELSE 

CALL  PRINTWd  . 'No'  .*500) 

END  IF 

10  CONTINUE 

END  IF 
CALL  CLRKB 
1  =  1 

20  CONTINUE 

IF(I  .EQ.  10)  I  =  1 

CALL  WLOCATEd  .IROW(  I  )  .ICOLd  )  .*500) 

IF(.NOT.  ISKEY(CHR.KEYTYPE) )  THEN 
GO  TO  20 
ELSE 

CALL  CTOA(CHR.IA.O) 

IF(KEYTYP  .EQ.  16  .AND.  lAd)  .EQ.  27)  GO  TO  30 
IF(KEYTYP  .EQ.  4  .AND.  lAd  )  .EQ.  13)  THEN 
1  =  1+1 
READ (*.1000)  CHR 
1000  FORMAT (Al) 

GO  TO  20 
ELSE 
END  IF 

;.F((KEYTYPE  .EQ.  2  .OR.  KEYTYPE  .EQ.  3)  .AND.  I  .EQ.  1)  THEN 
TITLE  =' 

WRITE(*.1200)  TITLE 

CALL  WLOCATEd  .IROWd  )  .ICOLd  )  .*500) 

READ(*.1100.ERR=20)  TITLE 

CALL  WLOCATEd  .IROWd  )  .ICOLd  )  ,*500) 

WRITE(*.1200)  TITLE 
1100  FORMAT (A40) 

12  00  FORMAT  dX.  A40,  $  ) 

1  =  1+1 
GO  TO  20 
ELSE 
END  IF 

IF(  (  (KEYTYPE.  EQ.  3  )  .AND.  ( IA(1 ) .  GE.  45.  AND.  lAd  ) .  LE.  57  )  )  .AND. 
1  (I  .EQ.  2  .OR.  I  .EQ.  3))  THEN 

CALL  PRINTWd  ,  '  '  .*500) 

CALL  WLOCATE(  1  .  IROW(  I  )  .ICOLd  )  .*500) 

READ(*.1300.ERR=20)  VAR 
1300  FORMAT(G10.  0) 

CALL  WLOCATEd  .IROWd  )  .ICOLd  )  .*500) 

IF(  VAR.  GE. -50000.  0  .AND.  VAR.  LE.  50000 .  0  )  VALUEd-l)  =  VAR 
WRITE(*.1400)  VALUEd-l) 

1400  FORMAT ( F8. 0 ,$ ) 

IF(VALUE(I-1)  .EQ.  VAR)  I  =•  I  +  1 
GO  TO  20 
ELSE 
END  IF 

IF(KEYTYP  .EQ.  2  .AND.  (I  .GE.  4  .AND.  I  .LE.  9))  THEN 
CALL  PRINTWd  .  '  '  .*500) 

CALL  WLOCATE(  1  ,  IROWd  )  .ICOLd  )  .*500  ) 

READ(*,1500.ERR=20)  IDUMMY 

A32 


n  n 


CALL  WLOCATE( 1 , IROW ( I ) . ICOL ( I ) . *500 ) 

1500  FORMAT (A3) 

IF(IDUMMY(1:1)  .EQ.  ‘Y*  .OR.  IDUMMY(1:1)  .EQ.  ‘y’)  THEN 
CALL  PRINTWd  ,  •  Yes-  .«500) 

FLAG (1-3)  =  .TRUE. 

1  =  1+1 
GO  TO  20 
£LS£ 

IFdDUMMYd  :1  )  .EQ.  'N'  .OR.  IDUMMYdd)  .EQ.  ’n-)  THEN 
CALL  PRINTWd.-No  ‘.»500) 

FLAG (1-3)  =  .FALSE. 

1  =  1+1 
GO  TO  20 
ELSE 

GO  TO  20 
END  IF 
END  IF 
ELSE 
END  IF 

CHR  =  GETCHX(KEYTYPE) 

CALL  WLOCATEd  .  IROW  ( I )  .  ICOL  { I ) .  »500  ) 

CALL  PRINTWd  .  '  '.*500) 

GO  TO  20 
END  IF 

30  CONTINUE 

CHR  =  GETCHX(KEYTYPE) 

CALL  CUROFF 

DIFFER  =  VALUE(2)  -  VALUE (1 ) 

IF((VALUE(1)  .LT.  -50000.  .OR.  VALUE(2)  .GT.  50000.)  .OR. 

IDIFFER  .LT.  1500.0)  THEN 
VALUEd)  =  -50000.0 
VALUE(2)  *  50000.0 
ELSE 
END  IF 

FMIN  =  VALUEd) 

FMAX  =  VALUE(2) 

500  RETURN 
END 

C 

c 

C  SOUND  PRESSURE  LEVEL  PLOTTING  PROGRAM 

C 

C 

SUBROUTINE  PLTSPL 

DIMENSION  XAR.;AY(203  )  ,  YARRAY(203) 

INTEGER*2  lOPORT . MODEL , NC . NPTS . NBLANK 
LOGICAL  FLAG(6) 

CHARACTER  CTEXT 
CWARACTER*40  TITLE 

COMMON  /RESULT/FLDNMR(201 ) .FLDNM(201 ) .FLEQ(201 ) .RANK(20.201 ) 
COMMON  /PLT /TITLE . FMIN . FMAX . FLAG . lOPORT . MODEL . FACT . I COUNT 

SET  PLOTTING  PARAMETERS 
FMIN  »  FLOAT ( IFIX ( FMIN/500. 0 ) )*500.0 
FMAX  =  FLOAT(IFIX(FMAX/500.0))»500.0 
DIFF  =  FMAX  -  FMIN 
DELTA  =  DIFF/ (FLOAT (10 ) ) 

AXLEN  =  -10.0 

IFMIN  »  IFIX(FMIN/500. 0)  +  101 
IFMAX  *  IFIX(FMAX/500.0)  +  101 
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NPTS  =  (IFMAX  -  IFMIN)  +  1 


IF(FLAG(1)'  THEN 
1  =  1 
ELSE 

IF(FLAG(2))  THEN 
1*2 
ELSE 

1  =  3 
END  IF 
END  IF 

0  CONTINUE 

CALL  PLOTS ( 0 . lOPORT . MODEL ) 

CALL  FONT(2) 

SET  FACTOR  FOR  HARD  COPY  DEVICE  OR  TERMINAL 
CALL  FACTOR(0.7) 

ONE  INCH  EQUALS  ONE  CENTIMETER  ON  HARD  COPY  DEVICE 
IFdCOUNT  .LE.  2)  CALL  FACTOR(  0 . 787  ) 


OTHER  DEVICE  SCALING  FACTOR 
IFdCOUNT  .EQ.  6)  CALL  FACTOR  (FACT) 

DRAW  THE  TITLE  AND  CENTER  IT  PLEASE 
NC  =  NBLANKC TITLE) 

X  =  6.0  -  (FLOAT(NC)/2.0)»(7. 5/40.0) 

CALL  SYMBOL(X.6.5,0.24.TITLE.0.0.40) 

DRAW  THE  AXIS  BOARDER 

CALL  STAXIS(-.177.  0.177,  .089,  .089,-1) 
CALL  AXISdl.  0.1. 0,  ■  ‘,0.5.0,90.0,0.0,0.0) 
CALL  AXISd.O.e.O,  '  •  .O.AXLEN, 0.0,0. 0,0. 0) 


SET  THE  DEFAULTS  TO  DRAW  ANNOTATED  AXIS 
CALL  STAXIS(  .177.  .177.  .089,  .089,  -1) 

DRAW  THE  X  AXIS 

CALL  AXIS (1.0, 1.0. 'Distance  in  f eet' . -16 , AXLEN , 0. 0 .FMIN .DELTA ) 


20 


NOW  USE  THE  FLAGS  TO  DRAW  THE  Y  AXIS 
IF(FLAG(5))  THEN 


CALL  AXISd. 0,1.0. 'Noise  Level '  ,  11  , -5 . 0 , 90 . 0 , 40 . 0 , 10 . 0  ) 


IF(FLAG(4))  THEN 
CTEXT  »  CHARd) 
IF(FLAG(1))  CALL 
IF(FLAGd))  CALL 
CTEXT  =  CHAR(2) 
IF(FLAG(2))  CALL 
IF(FLAG(2))  CALL 
CTEXT  =  CHAR'S) 
IF(FLAG(3))  CALL 
IF(FLAG(3))  CALL 
ELSE 
END  IF 
ELSE 

GO  TO  (20.  30.  40)  , 


SYMBOL( 9. 0, 5. 43 , 0. 15 .CTEXT ,0.0.-!) 
SYMBOL(9.5,5.33.0.15,  'Ldninr'  .0.0,5) 

SYMBOL(9.0.5.  10 . 0. 15 .CTEXT , 0. 0 .-1 ) 
SYMBOL(9.5.5. 00,0. 15. ‘ Ldn ‘ .0.0,3) 

SYMBOL(9.0,4.76.0. 15 . CTEXT , 0 . 0 , -1 ) 
SYMBOL (9. 5. 4. 66. 0.15. ‘ Leq ‘ .0.0,3) 


I 


CALL  AXISd. 0, 1. 0, ‘Ldnmr  in  dB'  .  11  ,-5. 0 . 90. 0 .40. 0 , 10. 0  ) 
CTEXT  *  CHARd) 
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IF(FLAG(4))  CALL  SYMBOL ( 9 . 0 . 5 . 43 . 0 . 15 . CTEXT . 0 . 0 . -1 ) 
IF(FLAG(4))  CALL  SYMBOL(  9 . 5 . 5 . 33 , 0.  15  .  '  Ldntnr  ‘  .  0.  0 . 5  ) 
GO  TO  50 

30  CALL  AXISd.O.l.O, ‘Ldn  in  dB'  .  9  . -5 .  0 . 90 .  0 , 40 . 0 . 10 .  0  ) 

CTEXT  =  CHAR(2) 

IF(FLAG(4))  CALL  SYMBOL ( 9 . 0 , 5 . 10 . 0 . 15 , CTEXT , 0 . 0 . -1 ) 
IF(FLAG(4))  CALL  SYMBOL( 9. 5 . 5. 00 . 0. 15 . ' Ldn '  . 0. 0 . 3  ) 

GO  TO  50 

40  CALL  AXISd.O.l.O, ‘Leq  in  dB' .  9 , -5 .  0 . 90 .  0 . 40 .  0 , 10.  0  ) 

CTEXT  =  CHAR(5) 

IF(FLAG(4))  CALL  SYMBOL( 9. 0 . 4 . 76 . 0. 15 , CTEXT , 0 . 0 , -1 ) 
IF(FLAG(4))  CALL  SYMBOL( 9. 5 , 4 . 66 . 0 , 15 , ' Leq '  , 0 . 0 , 3  ) 

END  IF 

50  CONTINUE 

C 

C  COPY  SOUND  PRESSURE  LEVEL  INTO  PLOTTING  ARRAY 

CALL  PLOTd.  0,1. 0,-3) 

60  K  =  0 

SUM  =  FMIN 

DO  100  J  =  IFMIN,  IFMAX 

K  *  K  +  1 

XARRAY(K)  =  SUM 
SUM  »  SUM  +  500 
GO  TO  (70.  80.  90) .  I 


70 

YARRAY ( K ) 
INTEO  =  1 
GO  TO  100 

»  FLDNMR(J) 

80 

YARRAY ( K ) 
INTEQ  *  2 
GO  TO  100 

*  FLDNM(J) 

90 

YARRAY ( K ) 
INTEQ  «  5 

*  FLEQ(J) 

100 

CONTINUE 

XARRAY(NPTS+1 )  -  FMIN 
XARRAY(NFTS+2 )  *  DELTA 
YARRAY(NPTS+1 )  *  40.0 
YARRAY(NPTS+2 )  *  10.0 

CALL  DRAW ( X ARRAY , YARRAY . NPTS , FLAG ( 4 ) . INTEQ ) 
IF(FLAG(5))  THEN 
1*1+1 

IF(FLAG(I)  .AND.  I  .LE.  3 )  GO  TO  60 
1*1+1 

IF(FLAG(I)  .AND.  I  .LE.  3 )  GO  TO  60 
CALL  PLOT (0. 0,0. 0.999) 

ELSE 
END  IF 

IF(.NOT.  FLAG(5))  THEN 

CALL  plot; 0.0, 0.0, 999) 

I  =  J  +  1 

IFCFLAGd)  .AND.  I  .LE.  3)  GO  TO  10 
1*1+1 

IF(FLAG(I)  .AND.  I  .LE.  3 )  GO  TO  10 

ELSE 
END  IF 
RETURN 
END 
C 
C 

C  PROBABILITY  OF  HIGHLY  ANNOYED  PLOTTING  PROGRAM 
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SUBROUTINE  PLTAPP 

DIMENSION  XARRAY(203).  YARRAY(203) 

INTEGER»2  lOPORT . MODEL . NC . NPTS . NBLANK 
LOGICAL  FLAG(6) 

CHARACTER  CTEXT 
CHARACTER»f40  TITLE 

COMMON  /RESULT/FLDNMR(201 ) .FLDNM(201 ) ,FLEQ(201 ) ,RANK(20 ,201 ) 
COMMON  /PLT /TITLE . FMIN , FMAX . FLAG . lOPORT . MODEL , FACT , I COUNT 
C 

C  SET  PLOTTING  PARAMETERS 

CALL  CLS 

FMIN  =  FLOAT(IFIX(FMIN/500. 0) )»500. 0 
FMAX  =  FLOAT ( IFIX ( FMAX/500. 0 ) )«500. 0 
DIFF  =  FMAX  -  FMIN 
DELTA  =  DIFF/ ( FLOAT ( 10 ) ) 

AXLEN  =  -10.0 

IFMIN  =  IFIX(FMIN/500. 0)  +  101 
IFMAX  =  IFIX ( FMAX/500. 0)  +101 
NPTS  =  (IFMAX  -  IFMIN)  +  1 
C 
C 

CALL  PLOTS ( 0 , lOPORT . MODEL ) 

CALL  FONT(2) 

C 

C  SET  FACTOR  FOR  HARD  COPY  DEVICE  OR  TERMINAL 
CALL  FACTOR(0.7) 

C 

C  ONE  INCH  EQUALS  ONE  CENTIMETER  ON  HARD  COPY  DEVICE 
IFdCOUNT  .LE.  2)  CALL  FACTOR(  0.  787 ) 

C 

C  OTHER  DEVICE  SCALING  FACTOR 

IFdCOUNT  .EQ.  6)  CALL  FACTOR  (FACT) 

C 

C  DRAW  THE  TITLE  AND  CENTER  IT  PLEASE 
NC  =  NBLANK (TITLE) 

X  *  6.0  -  (FLOAT(NC)/2.0)«(7. 5/40.0) 

CALL  SYMBOL(X ,6. 5,0.24 .TITLE, 0.0,40) 

C 

C  DRAW  THE  AXIS  BOARDER 

CALL  STAXIS(-. 177,  0.177,  .089,  .089,-1) 

CALL  AXISdi.O.l.O,  ■  ■.0,5.0.90.0.0.0.0.0) 

CALL  AXISd.  0,6.  0.  '  ',  0  .AXLEN ,  0.  0 , 0.  0 , 0 . 0 ) 

C 

C  SET  THE  DEFAULTS  TO  DRAW  ANNOTATED  AXIS 
CALL  STAXIS(  .177,  .177,  .089,  .089,  -1) 

C 

C  DRAW  THE  X  AXIS 

CALL  AXISd. 0,1.0, 'Distance  in  f  eef  , -16 ,  AXLEN ,  0. 0  .FMIN .  DELTA ) 
C 

C  COPY  SOUND  PRESSURE  LEVEL  INTO  PLOTTING  ARRAY 
K  »  0 

SUM  »  FMIN 
VALUE  *0.0 

DO  10  J  «  IFMIN,  IFMAX 
K  ■  K  +  1 
X ARRAY (K)  »  SUM 
SUM  «  SUM  +  500 

YARRAY(K)  -  1 . 0/(1 , 0+EXP( 10. 43-0. 132»FLDNMR(J )) ) 
IF(YARRAY(K)  .GT.  VALUE)  VALUE  «  YARRAY(K) 
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10  CONTINUE 

XARRAY(NPTS+1)  =  FMIN 
XARRAY(NFTS+2 )  =  DELTA 
YARRAY(NPTS+1 )  =0.0 

SET  THE  DEFAULTS  TO  DRAW  ANNOTATED  AXIS 
CALL  STAXIS(  .177,  .177,  .089,  .089,  1) 
IF( VALUE  .LT.  0.5)  THEN 
YARRAY(NPTS+2 )  *  0.1 


DRAW  THE  Y  AXIS 

CALL  AXIS ( 1 . 0 , 1 . 0 , ■ Prob.  of  High  Annoyance ’, 23 , -5 . 0 , 90 . 0 , 0 . 0 , 0 . 1 ) 
ELSE 

YARRAY(NPTS+2 )  =  0.2 
DRAW  THE  Y  AXIS 

CALL  AXISd. 0,1.0, 'Prob.  of  High  Annoyance '  ,  23  , -5 . 0 , 90 . 0 , 0 . 0 , 0 . 2  ) 
END  IF 

CALL  PLOT (1.0, 1.0, -3) 

CALL  DRAW ( X ARRAY , Y ARRAY , NPTS , . FALSE. , 1 ) 

CALL  PLOT (0. 0,0. 0,999) 

RETURN 

END 
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SUBROUTINE  RANKODR 
SINCLUDE:  ' spinattr . hdr ‘ 

INTEGER*2  KEYTYPE 
INTEGER  IA(1) 

LOGICAL  ISKEY,  FLAG (20) 

CHARACTER*!  CHR,  GETCHX ,  IDUMMY 
CHARACTER*2  OPCR(300) 

CHARACTER*3  ACC(300),  AIRSPD(300) 

CHARACTER'S  POWSET(300) 

CHARACTER*6  POWUNTOOO) 

CHARACTER*?  NA!V!E(300) 

CHARACTER*20  POWDES(300) 

CHARACTER*25  DRAG (300) 

CHARACTER*50  OUTPUT 
CHARACTER*55  TEXT(6) 

COMMON  /AIRPARA/ALT(20) .  FMEAN(20).  ICODE(20).  IREL.  I TEMP . 

1  MONTH,  NUMAIR.  FNUMD(20).  FNUMN(20).  POWER(20).  SEL(20.25). 

2  SIGMA,  VEL(20) 

COMMON  /FRAME/SCN (1024) 

COMMON  /NOISE/MAX . ACC , AIRSPD . DRAG , NAME , OPCR , POWDES , POWSET , POWUNT 
COMMON  /RESULT/FLDNMR(201 ) .FLDNM(201 ) .FLEQ(201 ) ,RANK(20,201 ) 

DATA  TEXT(l)  /'Enter  the  distance  in  feet  from  corridor  centerline 
1  •  / 

DATA  TEXT(2)  /'to  ground  location  point  of  interest  :'/ 

DATA  TEXT ( 3 )  /'Distance  Out  Of  Range'/ 

DATA  TEXT(4)  /'  RanX  Contributors  At  Distance  ®'/ 

DATA  TEXT(5)  /'  Config  Aircraft  Ldnmr  Prob  of  High 

lly'  / 

DATA  TEXT(6)  /'  Number  Name  Annoyed  ' 

1/ 

C 

DO  10  I  «  1.  20 

FLAG(I)  =  .TRUE. 

10  CONTINUE 

20  CALL  DEFWINDd  .8,0,14,79,DBLBDR.NORMAL) 

CALL  WINDBUF(SCN) 

CALL  CLS 

CALL  CAPTSCN(SCN) 

CALL  CLRB(l) 

CALL  BLDBd  ,1 ,10, TEXTd  )  .0**500) 

CALL  BLDBd  ,2 .10.TEXT(2  )  .0.*500) 

CALL  DISPSCN(SCN) 

CALL  CURON 
GO  TO  40 

30  READ(*,1000)  IDUMMY 

1000  FORMAT (Al) 

CALL  WLOCATE(1.2.48.*500) 

CALL  PRINTWd  .  '  ',*500) 

40  CALL  WLOCATEd  .2 ,48.*500) 

IF (.NOT.  ISKEY (CHR. KEYTYPE ) )  THEN 
GO  TO  40 
ELSE 

CALL  CTOA(CHR.IA.O) 

IF( (KEYTYPE  .EQ.  3)  .AND. 

1  ((IA(1).GE.45  .AND.  IA(1).LE.57)  .AND.  IA(1).NE.47))  THEN 

READ(*.1100,ERR»30)  DIST 
WRITE(OUTPUT,1100) 

1100  FORMAT(F10. 0) 

ELSE 

CHR  *  GETCHX ( KEYTYPE ) 
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CALL  WLOCATEd  .2 ,48.»500) 

CALL  PRINTWd,-  *.*500) 

GO  TO  40 
END  IF 
END  IF 
CALL  CUROFF 

IF(DIST  .LT.  -50000.0  .OR.  DIST  .GT.  50000.0)  THEN 
CALL  DEFWINDd  ,8,0  d4. 79.  DBLBDR,  NORMAL) 

CALL  WINDBUF(SCN) 

CALL  CLRBd) 

CALL  BLDBd  ,2  .20,TEXT(3)  ,0.»500) 

CALL  DISPSCN(SCN) 

CALL  PAUSE 
GO  TO  20 
ELSE 
END  IF 

DETERMINE  WHICH  VECTOR  WAS  PICKED 
DIST  *=  DIST  -t-  50000.0 
IVAL  =  IFIX(DIST/500.0)  +  1 
P  =  MOD(DIST/500. 0.1. 0) 

IF(P  .GT.  0.5)  IVAL  »  IVAL  +  1 

PREPARE  SCREEN 

CALL  DEFWIND (2,0.11,24.68. DBLBDR . NORMAL ) 

CALL  WINDBUF(SCN) 

CALL  CLS 

CALL  CAPTSCN(SCN) 

CALL  CLRB(2) 

CALL  BLDB(2 .0.1 .TEXT(4) ,0.»500) 

CALL  BLDB(2 .1 .1 .TEXT(5) .0.»500) 

CALL  BLDB(2 .2 .1 ,TEXT(6) .0.»500) 

CALL  DISPSCN(SCN) 

DIST  =  (IVAL  -  101)*500.0 
WRITE(OUTPUT,1200)  DIST 
FORMAT ( G8. 3E1 ) 

CALL  WLOCATE(2,0,40,»500) 

CALL  PRINTW(2  .OUTPUTd  :8)  .»500) 

DO  80  I  *  1.  NUMAIR 
J  *  0 
J  *  J  +  1 

IF(.NOT.  FLAG(J)  .AND.  J  .LT.  NUMAIR)  GO  TO  50 
IF(J  .EQ.  NUMAIR)  THEN 
ITRACK  =  NUMAIR 
GO  TO  70 
ELSE 
END  IF 

VALMAX  =  RANK (J, IVAL) 

ITRACK  =  J 

DO  60  K  =  1.  NUMAIR 

IF( RANK (K, IVAL)  .GT.  VALMAX  .AND.  FLAG(K))  THEN 
VALMAX  *  RANK (K, IVAL) 

ITRACK  =  K 
ELSE 
END  IF 
CONTINUE 

FLAG (ITRACK)  =  .FALSE. 

ANNOY  =  1.0/ (1.0  +  EXP (10. 43  -  0. 132*RANK( ITRACK . IVAL )) ) 
CALL  WLOCATE(2 .1+2 ,2 ,*500) 

WRITE(OUTPUT.1300)  ITRACK . NAME ( I CODE ( ITRACK ) ) . 
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1  I FIX (RANK (I TRACK. I V AL) ) .ANNOY 

1300  FORMAT ( 4X . 1 3 . 9X . A7 . 5X . 1 3 , 9X . F6 . 3 ) 
CALL  PRINTW(2  .OUTPUTd  :50)  .*500) 
80  CONTINUE 

CALL  ANYKEY 
500  RETURN 
END 
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c 

c 

C  ADDS  TWO  NUMBEES  IN  DECIBELS 

C 

C 

SUBROUTINE  ADDER 
SINCLUDE;  ’ spinattr , hdr ' 

DIMENSION  IA(1).  ICOL(3).  VALUEO) 

INTEGER*2  KEYTYPE 
LOGICAL  ISKEY 

CHARACTER  IDUMMY.  CHR .  GETCHX 
CHARACTER*34  TEXT(5) 

COMMON  /FRAME/SCN(1024) 

PARAMETER  (IROW=6) 

DATA  TEXTd  )/■  DECIBEL  ADDITION'/ 

DATA  TEXT (2 )/' Enter  Noise  Levels  In  Decibels'/ 

DATA  TEXT(3)/'  dB  +  dB  *  dB'/ 

DATA  TEXT ( 4 )/' press  RETURN  to  advance  the  cursor'/ 

DATA  TEXT(5)/'press  ESC  when  all  done'/ 

DATA  ICOL/10.  21.  32/ 

DATA  VALUE/3«0.0/ 

CALL  CURON 

CALL  DEFWIND (1,5.15.17,65. DBLBDR . NORMAL ) 

CALL  WINDBUF(SCN) 

CALL  CLS 

CALL  CAPTSCN(SCN) 

CALL  CLRB(l) 

CALLBLDBd.  1 , 1 7  ,  TEXT  (1  )  .  0  .  »100 ) 

CALL  BLDBCl.  3 . 10 . TEXT ( 2 ) . 0 .•100 ) 

CALLBLDBd.  6 . 10  .  TEXT  (  3  )  .  0  .  »100  ) 

CALL  BLDBd.  9.  8  .TEXT<  4  )  .0  .«100) 

CALL  BLDBd  .10,13.TEXT(5)  .0.«100) 

CALL  DISPSCN(SCN) 

C 

C  INPUT  RESULTS,  ALLOW  FOR  EDITING  FEATURES 

1  =  1 
GO  TO  20 

10  READ (». 1000)  IDUMMY 

1000  FORMAT (Al) 

CALL  WLOCATEd  .IROW.ICOLd  )  .•100) 

CALL  PRINTWd  ,  '  '.•lOO) 

20  CONTINUE 

IF(I  .EQ.  3)  THEN 

CALL  WLOCATEd.IROW.ICOL(I),^100) 

1=1 

VALUE(3)  =  10. 0^^(VALUEd  )/10. 0)  +  10 . 0*^  ( VALUE(  2  ) /lO .  0  ) 
IF( VALUEO)  .LE.  0.0)  THEN 

CALL  PRINTWd  ,'••••*'  ,•100  ) 

ELSE 

VALUE(3)  =  10. 0^LOG10(VALUE(3 ) ) 

WRITE(^.llOO)  VALUEO) 

1100  FORMAT ( F5. 1 .S ) 

END  IF 
ELSE 
END  IF 

CALL  WLOCATEd  .IROW.ICOLd  )  ,•100) 

IF (.NOT.  ISKEY(CHR.KEYTYPE) )  THEN 
GO  TO  20 
ELSE 

CALL  CTOA(CHR.IA,0) 
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1200 

1300 

30 

100 


IF(KEYTYPE  .EQ.  16  .AND.  IA(1)  .EQ.  27)  GO  TO  30 
IF(KEYTYPE  .EQ.  4  .AND.  lA(l)  . EQ.  13)  THEN 
1  =  1+1 

CHS  =  GETCHX ( KEYTYPE ) 

GO  TO  20 
ELSE 
END  IF 

IF( (KEYTYPE  .EQ.  3)  .AND. 

1  ( (IA(1 ).GE.45  .AND.  IA(1).LE.57)  .AND.  IA( 1 ) . NE. 47 ) )  THEN 

CALL  PRINTWd  .  '  -,*100) 

CALL  WLOCATE( 1 . I ROW . I COL ( I ) . *100 ) 

READ (*, 1200. ERR=10)  VALUE(I) 

IF(VALUE(I)  .GT.  140.0  .OR.  VALUE(I)  .LT.  0.0)  GO  TO  20 
CALL  WLOCATE( 1 . I ROW , ICOL( I ) . *100 ) 

WRITE(*.1300)  VALUE(I) 

FORMATCFS.O) 

FORMAT ( F6 . 1  .  $ ) 

1  =  1+1 
GO  TO  20 
ELSE 

CHR  =  GETCHX (KEYTYPE) 

CALL  WLOCATEd  .IROW.ICOL(I)  .*100) 

CALL  PRINTWd.'  '.*100) 

GO  TO  20 
END  IF 
END  IF 
CONTINUE 

CHR  =  GETCHX (KEYTYPE) 

CALL  CUROFF 

RETURN 

END 
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This  is  the  Header  file  for  specifying  Screen  Attributes 
and  Window  Parameters 


The  Color  Attributes  are  of  2  types:  Those  determining  the 
BACKGROUND  Color,  and  those  determining  the  CHARACTER  Color. 

The  Colors  are  BlacK  (BLK),  Dark  Blue  ( DBLU ) .  Green  (GREEN), 
Red  (RED),  Lavender  (LAV).  Yellow  (YELLO).  and  White  (WHITE). 

To  make  a  particular  Color  combination  on  the  screen,  merely 
add  up  the  two  colors  desired  for  the  BACKGROUND  ( BKG ) ,  and 
the  CHARACTERS  (CHR). 

For  example,  to  place  RED  characters  on  a  GREEN  background, 
use  the  sum:  REDCHR  +  GREENBKG 

If  the  Characters  are  to  be  of  NORMAL  Intensity,  Nothing  else 
need  be  done.  If  the  Characters  are  to  be  BRITE  also,  use 
REDCHR  +  GREENBKG  +  BRITE 

If  they  are  to  FLASH  also,  use: 

REDCHR  +  GREENBKG  +  BRITE  +  FLASH 


The  WINDOW  Parameters  are  SCROLL,  WRAP,  DBLBDR  for  Double  Border, 
SINGED  for  Single  Border,  LIST  for  no  carriage  control, 
and  NONE  for  no  Parameters 

If  a  Window  was  to  have  Scrolling  and  a  Double  Border,  use 
SCROLL  +  DBLBDR 

Note  that  if  LIST  is  not  specified,  then  Fortran  will  assume 
that  normal  formatting  rules  are  in  effect:  The  first  character 
will  be  assumed  to  be  a  Carriage  Control  Character. 


INTEGER  BLKBKG,  BLKCHR ,  DBLUBKG.  DBLUCHR .  GREENBKG,  GREENCHR . 
LBLUBKG,  LBLUCHR,  REDBKG,  REDCHR.  LAVBKG .  LAVCHR . 

YELLOBKG,  YELLOCHR ,  WHITEBKG.  WHITECHR,  FLASH,  UNDRLINF, 
NORMAL,  BRITE,  BRIGHT.  INVERSE.  SCROLL.  WRAP,  DBLBDR,  SINGLBD, 
UNITZERO,  DEFAULT 


PARAMETER 

(BLKCHR  = 

0) 

PARAMETER 

(BLKBKG  = 

0) 

PARAMETER 

(DBLUBKG  = 

:  16) 

PARAMETER 

(DBLUCHR  = 

^  1  ) 

PARAMETER 

( GREENBKG 

=  32  ) 

PARAMETER 

(GREENCHR 

=  2) 

PARAMETER 

(REDBKG  = 

64  ) 

PARAMETER 

(REDCHR  = 

4) 

PARAMETER 

(LAVBKG  = 

80  ) 

PARAMETER 

(LAVCHR  = 

5) 

PARAMETER 

(YELLOBKG 

=  96) 

PARAMETER 

(YELLOCHR 

=  6) 

PARAMETER 

(WHITEBKG 

=  112  ) 

PARAMETER 

(WHITECHR 

=  7) 

PARAMETER 

(UNDRLINE 

=  1  ) 

PARAMETER 

(DEFAULT  = 

^  0) 

PARAMETER 

(NORMAL  = 

7  ) 
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PARAMETER  (FLASH  =  128) 
PARAMETER  (INVERSE  =  112) 
PARAMETER  (BRIGHT  =  15) 
PARAMETER  (BRITE  =  8) 
PARAMETER  (LBLUBKG  =  48) 
PARAMETER  (LBLUCHR  =  3) 
PARAMETER  (NONE  =  0) 
PARAMETER  (SINGLED  =  1) 
PARAMETER  (UNITZERO  =  0) 
PARAMETER  (DBLBDR  =  3) 
PARAMETER  ( WRAP  =  4  ) 
PARAMETER  (SCROLL  =  8) 
PARAMETER  (LIST  =  16) 
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